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ABSTRACT 

We show that it is possible to construct e'/e to NLO using partially 
quenched chiral perturbation theory (PQChPT) from amplitudes 
that are computable on the lattice. We demonstrate that none of 
the needed amplitudes require three-momentum on the lattice for ei- 
ther the full theory or the partially quenched theory; non-degenerate 
quark masses suffice. Furthermore, we find that the electro-weak 
penguin (A7 = 3/2 and 1/2) contributions to e'/e in PQChPT can 
be determined to NLO using only degenerate (m K = m„) K — > it 
computations without momentum insertion. Issues pertaining to 
power divergent contributions, originating from mixing with lower 
dimensional operators, are addressed. Direct calculations of K — > tttt 
at unphysical kinematics are plagued with enhanced finite volume ef- 
fects in the (partially) quenched theory, but in simulations when the 
sea quark mass is equal to the up and down quark mass the enhanced 
finite volume effects vanish to NLO in PQChPT. In embedding the 
QCD penguin left-right operator onto PQChPT an ambiguity arises, 
as first emphasized by Golterman and Pallante. With one version 
(the "PQS") of the QCD penguin, the inputs needed from the lattice 
for constructing K — > tttt at NLO in PQChPT coincide with those 
needed for the full theory. Explicit expressions for the finite loga- 
rithms emerging from our NLO analysis to the above amplitudes are 
also given. 



1 Introduction 

There have been several recent lattice attempts to calculate Re(e'/e), the direct 
CP violating parameter in K — ► tttt decays. These include attempts with do- 
main wall fermions by the CP-PACS [T] and RBC g] Collaborations. A notable 
feature of both of these calculations is that their central values differ drastically 
from experiment. The experiments at CERN 3 and Fermilab [3] have yielded 
an experimental grand average of Re(e'/e) = (1.8 ± 0.4) x 10~ 3 jS]. The lattice 
collaborations find a value ~ —0.5 x 10 -3 , a negative value, though the groups 
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have made rather severe approximations. Such a disagreement between theory 
and experiment should not be totally unexpected given the serious approxima- 
tions and resulting systematic errors, which have so far been necessary in order 
to implement the calculation on the lattice. |SJ[7] 

One of these uncontrolled approximations was the use of the quenched ap- 
proximation, where the fermion determinant in the path integral is set to a 
constant in order to make the problem more tractable on present day comput- 
ers. Another was the use of leading order chiral perturbation theory (ChPT) 
to relate unphysical K — ► tt and K — > amplitudes to the physical K — > tttt 
amplitudes, as first proposed by Because of the difficulty of extracting mul- 
tihadron decay amplitudes from the lattice, as expressed by the Maiani- Testa 
theorem UJ, it is much easier to compute the two- and three-point functions 
(i.e., K — > and K — > tt, respectively) and use ChPT to extrapolate to the 
physical matrix elements. 

It is likely that the next-to-leading order (NLO) corrections to ChPT will 
be significant for the operators that contribute to Re(e'/e), and should not be 
neglected. Unfortunately, at higher orders in ChPT the number of free pa- 
rameters that must be determined from first-principles methods like the lattice 
proliferates rapidly. It has been shown by Cirigliano and Golowich ^U] that the 
dominant electroweak penguin contributions [(8,8) 's] to K — > tttt can be recov- 
ered at NLO from K — ► tt amplitudes using 4-momentum insertion. Bijnens, et 
al. m showed how to obtain most of the low-energy constants (LEC's) relevant 
for the case of the (8,l)'s and (27,l)'s using off-shell K — ► tt Green's functions; 
not all LEC's could be determined using this method, though. 

In ^2], it was shown how to obtain physical K — * tttt, AI = 3/2 [(27,l)'s and 
(8,8) 's] at NLO from K — > tttt at unphysical (SPQcdR) kinematics accessible to 
the lattice. This method requires 3-momentum insertion, and it is not yet clear 
if it can be extended to the AI = 1/2, K — > tttt amplitudes. In our previous 
paper |13j . an alternative method was proposed for constructing the physical 
K — ► tttt amplitudes to NLO for all (AI =1/2 and 3/2) operators of interest. 
For the AI = 3/2 amplitudes this requires K — > K; K —* tt, AI = 3/2; and 
K — ► tttt, AI = 3/2 at one of (at least) two unphysical kinematics points where 
the Maiani- Testa theorem can be bypassed. The two special kinematics points 
where this is possible have been discussed in the literature: (i) rrijf = m^*, 
where the weak operator inserts energy |14| : and (ii) m 1 ^ = 2mJ r " i , i.e. at 
threshold [E]. As in J^jj we refer to these two cases as unphysical kinematics 
point 1 (UK1) and point 2 (UK2), respectively. Finally, it was also shown in 
how to obtain the physical K — > tttt at NLO for the AI = 1/2, (8,1) (e.g. 
Qn and Qe) and the mixed (27, 1) © (8, 1) case (e.g., Q2) using K — > tt with 
4-momentum insertion and K — ► tttt at both UK1 and UK2. Note that the 
mixed case also requires information obtainable from the amplitudes needed to 
get K — > tttt, AI — 3/2 for the (27,l)'s. The main purpose of ^Sj was, in fact, to 
show that even for the (8,l)'s all of the information needed to construct K ^ tttt 
to NLO in ChPT could be obtained from amplitudes that can be computed on 
the lattice, at least in principle. 

There are other unphysical kinematics values for the K — ► tttt amplitudes 
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where the initial and final state mesons are at rest that bypass the Maiani- Testa 
theorem. These kinematics are similar to UK1 in that they require energy inser- 
tion, but with itik ^ m^. We call this set of kinematics UKX. This corresponds 
to the SPQcdR kinematics with both pions at rest (T2|. Lattice calculations at 
these values of the kinematics are important given that the calculation at UK1 
has difficulties and also given that it will be important to determine the 
NLO LEC's in as many ways as possible for additional redundancy. Even if it 
is difficult or impossible to obtain the necessary NLO low energy constants for 
the (8,l)'s at UK1, one can obtain the same information using UKX. Thus, all 
information for the (8,l)'s can be determined to NLO without using UK1, the 
difficulties of which are discussed in Section 8 and in the note added in revision. 
Results at UKX are also given in Section 8. 

In this work we show that where 4-momcntum insertion is required for any 
of the amplitudes needed according to the prescription of it suffices to allow 
only energy insertion at the weak operator such that the initial and final state 
mesons are at rest. This means that the K — > tttt amplitudes can be constructed 
to NLO using non-degenerate quarks, but without using 3-momentum insertion, 
making the computation much more economical. 

Another approach to K — * tttt and e'/e amplitudes has been proposed by 
Lcllouch and Luscher |17| in which finite volume correlation functions on the 
lattice are used to extract physical amplitudes without recourse to ChPT, at 
least in principle. This method is expected to be difficult computationally, but a 
way of reducing the cost of the Lellouch-Luscher method has been proposed |T8). 
An alternative method to obtain K irir amplitudes to all orders in ChPT has 
been proposed by ^Hj; this proposal makes use of dispersion relations. Both of 
the above methods depend crucially on unitarity, so it is unclear if they can be 
implemented with partially quenched lattice simulations. 

Although NLO ChPT may not be the final answer, it is more reliable than 
leading order, and it is useful to have the NLO expressions even to extract the 
leading order LEC's from the lattice data. Since the lattice data that will be 
generated in the near term will be in the (partially) quenched approximation, it 
is necessary to have the corresponding amplitudes in partially quenched ChPT. 
Therefore, in this paper, we present the partially quenched expressions for the 
quantities of greatest interest for Re(e'/e), namely the amplitudes for the (8,1) 
and (8,8) operators. For the partially quenched amplitudes we assume that all 
relevant quark masses are small compared to the mass, so that the r/' can be 
integrated out, and the LEC's of the partially quenched theory coincide with 
those of the full theory when the number of sea quarks is three |2(J| . 

For the AI —1/2 amplitudes there is an additional complication involving 
eye diagrams having to do with the sum over quarks in the penguin operators 
|21|. For the left-right gluonic penguin operators the two possible choices cor- 
respond to what we will call the PQS (partially quenched singlet) method and 
the PQN (partially quenched non-singlet) method. They are discussed in detail 
in Section 6.1. It is important to note that only for the PQS method can the 
LEC's sufficient to construct e'/e to NLO be determined, whereas it is not clear 
if the PQN method can be extended to NLO. Indeed, a significant advantage 
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of the PQS implementation is that the ingredients needed from the lattice to 
obtain all K — > tttt amplitudes to NLO in PQChPT are the same as in the full 
theory. Therefore, the PQS method is used to compute the NLO amplitudes in 
this paper. Finally, it should be mentioned that the AI = 1/2, K — > tttt am- 
plitudes receive enhanced finite volume contributions in the partially quenched 
theory |221 I23| . However, when m sea — m u = rrid, the infra-red divergences in 
the K — > tttt amplitudes (at UK1 and UK2) vanish in PQChPT in the infinite 
volume Minkowski space amplitudes. In an earlier version of this paper we had 
pointed out that it would be important to study the finite volume effects of 
these amplitudes; the corresponding finite volume Euclidian Green's functions 
were calculated by |15| while this work was in revision. 1 

In the partially quenched theory, it is possible to construct the (8,8) K — > tttt 
amplitudes to NLO using only degenerate valence quark masses in K — > tt, along 
with K — ► in order to perform the power divergent subtraction in the AI =1/2 
case. Additional redundancy is possible if one uses nondegenerate valence quark 
masses in the K — > tt calculation. 

The content of the paper is as follows. Section 2 briefly reviews the formalism 
of effective four-fermion operators in a standard model calculation. Section 3 
reviews ChPT and the realization of the effective four-quark operators in terms 
of ChPT operators for weak processes. Section 4 reviews partially quenched 
chiral perturbation theory and how it can be extended to the electroweak sector. 
Section 5 presents results for the full theory, demonstrating that for all the 
amplitudes considered in |13j. 3- momentum insertion is not essential and non- 
degenerate quark masses suffices to construct K — > tttt to NLO. In Section 6 
a discussion of the treatment of eye-diagrams in the partially quenched theory 
is given, as well as a comparison of PQS and PQN results at leading order 
according to the papers by Golterman and Pallante [2] • Sections 7 and 8 present 
the main results of this paper, showing how to obtain the K — > tttt amplitudes 
needed for Re(e'/e) in the partially quenched theory from quantities which can 
be computed directly on the lattice. Section 7 deals with the (8,8) amplitudes, 
while Section 8 deals with the (8,l)'s. Section 9 discusses the checks done on 
the various one-loop logarithmic expressions. Section 10 presents the conclusion. 
Section 11 is a note added in revision. The finite logarithm contributions to the 
relevant amplitudes are presented in a set of Appendixes. Errors in Eqs (31, 
D6) of ^Sj are corrected in Appendix F. 

2 Effective Four Quark Operators 

In the Standard Model, the nonleptonic interactions can be expressed in terms of 
an effective AS = 1 hamiltonian using the operator product expansion |2"4"ll2"K] , 

Ref 1151 found that the infra-red problems do not vanish for UK1 finite volume Euclidean 
correlation functions in the partially quenched theory; for further details, see our note added 
in revision. 
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(HWas=i|^) = ^J2 V CKMCi(^(nn\Q i \K) fl , (1) 

where Vq km are the relevant combinations of CKM matrix elements, Ci(/i) are 
the Wilson coefficients containing the short distance perturbative physics, and 
the matrix elements I^k^Q^K) ^ must be calculated nonperturbatively. The 
four quark operators are 



'l = S a Jfj,(l 


-7 5 )^aU 6 7 Al (l-7 5 )wf>, 


(2) 


] 2 = SoTmC 1 


-7 5 )4%7"(1 ~7 5 K, 


(3) 


= 807,4(1 ~ 


7 5 )d Q ]Tg b y<(l- 7 5 )q b , 


(4) 


= s 7 M (l - 


T 5 )4^| T "(l-7 5 ) ? „, 


(5) 




7 5 )d Q ^g b7 ^(l+7 5 )%, 
9 


(6) 


= S a 7p(l ~ 


7 5 HE^(l+7 5 )?a, 
9 


(7) 




7 5 R J] e g g b 7 p (l + 7 5 )%, 
9 


(8) 


3 - n 


7 5 H^e g g b7 ^(l + 7 5 )g Q , 
9 


(9) 


2Sa7^(!- 


7 5 K^e 9 g b7 Al (l-7 5 )'Z6, 
9 


(10) 




-7 5 H^e g g b7 ^(l- 7 5 )g Q . 


(11) 



In the effective theory Qi and Q2 are the current-current weak operators, 
Q3 ~ Qe are the operators arising from QCD penguin diagrams, while Q7 — Q10 
are the operators arising from electroweak penguin diagrams. Note that the 
definitions of Q\ and Q2 are different from our previous paper |13|. After a 
Fierz transformation, one can see that the definitions of the two operators are 
switched. We have changed the definitions of Q\ and Q2 to be consistent with 
the basis used by RBC [5J and that of [21]; this does not, of course, effect any 
of the results of our previous paper. 
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3 Chiral Perturbation Theory 



Chiral perturbation theory (ChPT) is an effective quantum field theory where 
the quark and gluon degrees of freedom have been integrated out, and is ex- 
pressed only in terms of the lowest mass pseudoscalar mesons [2S| - It is a 
perturbative expansion about small quark masses and small momentum of the 
low mass pseudoscalars. The effective Lagrangian is made up of complicated 
nonlinear functions of the pseudoscalar fields, and is nonrenormalizable, making 
it necessary to introduce arbitrary constants at each order in perturbation the- 
ory. In such an expansion, operators of higher order in the momentum (terms 
with increasing numbers of derivatives) or mass appear at higher order in the 
perturbative expansion. The most general set of operators at a given order can 
be constructed out of the unitary chiral matrix field S, given by 



exp 



2i(f) a \ a 



f 



(12) 



where A a are proportional to the Gell-Mann matrices with tr(A Q Ab) = 6 a b, <t> a 
are the real pseudoscalar-meson fields, and / is the meson decay constant in the 
chiral limit, with / ff equal to 130 MeV in our convention. 

At leading order [0(p 2 )] in ChPT, the strong Lagrangian is given by 

4? = ^tr[0 M E0"£] + -^-tr^E + E t x]j (13) 
where x = ( m u, m d , m s ) d - mg and 

D = m l+ _ "4 + _ m K0 

m u -\-md m u -\-m s md+m s 

The leading order weak chiral Lagrangian is given by [Sj I10| 



d§ = a 88 tr[A 6 EQEt] + aitr[A 6 ^S^St] + a2 2B tr[A 6 ( x t S + E+x)] 

+a 27 t%(Xd^)UZd^) l 3 + H.c, (14) 

where is symmetric in i,j and k, I, traceless on any pair of upper and lower 
indices with nonzero elements = 1, ^22 = 1/2 and = —3/2. Also, Q is the 
quark charge matrix, Q = 1/3(2, — 1, — l)diag and {Xe)ij = <5i3^j2- The reason 
A6 enters these expressions is because it picks out the s to d, AS — 1 transition. 

The terms in the weak Lagrangian can be classified according to their chi- 
ral transformation properties under SU(3)l x SU(3)ij. The first term in (14) 
transforms as 8^ x 8r under chiral rotations and corresponds to the electroweak 
penguin operators Q7 and The next two terms in (14) transform as 8l X 1_r, 
while the last transforms as 27^ x 1^ under chiral rotations. All ten of the four 
quark operators of the effective weak Lagrangian have a realization in the chiral 
Lagrangian differing only in their transformation properties and the values of 
the low energy constants which contain the non-perturbative dynamics of the 
theory. 



G 



For the transition of interest, K — ► irir, the operators can induce a change 
in isospin of \ or | leading to a final isospin state of the pions of or 2, respec- 
tively. We can then classify the isospin components of the four quark operators 
according to their transformation properties ^ |2] : 





Q2 ■• Q§ -i ^ 


->l/2 . 
tflO • 


8 L x \ R @ 27 L x 1 R 


,3/2 


^2 1 ^9 i 1 


-,3/2 . 


27 L x 








Qe /2 : 8 L x 




£?7^ 1 Qs^ ' 


Q 3 / 2 , 


Qs /2 : 8l x 8 fl . 



Note that Q3 — <?6 are pure isospin | operators. At NLO the strong La- 
grangian involves 12 additional operators with undetermined coefficients. These 
were introduced by Gasser and Leutwyler in |27|. The complete basis of coun- 
terterm operators for the weak interactions with AS — 1,2 was treated by 
Kambor, Missimer and Wyler in [2E] and [21] • A minimal set of counterterm 
operators contributing to K — > n and K — » 7T7t for the (8l, 1/j) and (27l, 1_r) 
cases is given by jHO], with the effective Lagrangian 

4™°> ^ e ,of +$>of 7 ' 1) +E c A (M) > ( 15 ) 
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(8,1) 
(8,1) 
(8,1) 
(8,D 
(8,1) 

w ll 

/*)(8,1) 
1 - / 12 
/n(8,l) 
^13 

^15 

/n(8,l) 

^35 

^39 



0(8,8) 


= tr 


0(8,8) 


= tr 


0(8,8) 


= tr 


0(8,8) 


= tr 


0(8,8) 


= tr 


ro (8,8) 


= tr 



X 6 S]tr[S], 

A 6 P 2 ], 

A 6 P]tr[P], 

X 6 [S,P}}, 

A 6 {S,L 2 }], 

A 6 L M ]tr[{ZAS}], 
X 6 S][L% 
X 6 [P,L% 
\ & {L^d v W^}], 

XeL^QELv], 
A 6 L M ]tr[StQSL^], 
A 6 {Etg S ,L 2 }], 
Aa{EtQE,5}], 

A 6 [stgs,p]], 

A 6 EtgE]tr[S'], 



(27,1) 



(27,1) 



(27 
20 
(27, 



(27,1) _ 



(27,1) _ 



( 



tt?i(L^([i^,P])^ 



t^(^)*(^)Jtr[S], 



(16) 



with 5 = 2P (x f E + E+x), P = 2P (x f E - E+x), ^ = *E^<9 M E , and W^" = 
2{d li L v +d v L v ). 

This list is identical to that of Bijnens et al. for the (27, l)'s and the 
(8,l)'s, except for the inclusion of C^'gg and 0^2i which contain surface 
terms, and so cannot be absorbed into the other constants for processes which do 
not conserve 4-momentum at the weak vertex. Since we must use 4-momentum 
insertion in a number of our amplitudes, these counterterms must be considered, 
and they are left explicit even in the physical amplitudes. The list of (8,8) 
operators is that of Cirigliano and Golowich [TUJ . 

The divergences associated with the counterterms have been obtained in |1U|. 
|11| . and [2HI- The subtraction procedure can be defined as 



1 



167T 2 / 2 



' +^(7E-l-ln47r) 



d-4 2 



1 



dl ~ + 16rr 2 / 2 
c-e 1 

— Oj 



! + -(7b - 1 -ln47r) 



d-4 2 



16tt 2 / 2 



+1( 7B _ 1-1*4*0 



2(aie< + a^-), 

2«277i; 

2a 8 s?? l , 



(17) 
(18) 
(19) 
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with the divergent pieces, £i,e^,7i, f]% given in Table 1. 

It is also necessary for the method of this paper to consider the 0(p 4 ) strong 
Lagrangian, which was first given by Gasser and Leutwyler, = Y^LiO^. 



Table 1: The divergences in the weak 0(p 4 ) counterterms, e<'s and d^s, for 
the (8,l)'s and (27,l)'s, respectively, and the divergences in the weak 0(p 2 ) 
counterterms, the Cj's for the (8,8) 's. 





e* 


< 


di 


7i 


Ci 




1 


1/4 


5/6 


1 


-1/6 


1 





2 


-13/18 


11/18 


2 





2 


-2 


3 


5/12 





4 


3 


3 


-3/2 


4 


-5/36 





5 


1 


4 


3/2 


5 





5/12 


6 


-3/2 


5 





10 


19/24 


3/4 


7 


1 


6 


1 


11 


3/4 





20 


1/2 






12 


1/8 





24 


1/8 






13 


-7/8 


1/2 










15 


23/24 


-3/4 










35 


-3/8 













39 


-3/16 














The strong 0(p 4 ) operators relevant for this calculation are the following 27 : 



Q {st) 


-tr[Z/f, 


o ( 2 st) 


= tr[L li L„]tr[L»L v 


(st) 


= tr[L 2 L 2 }, 


Q (st) 


= tr[L 2 ]tr[S'], 


Q {st) 


= tr[L 2 S], 


n(st) 

^6 


- tr[5] 2 , 


n {st) 
^8 


= \tr[S 2 -P 2 ]. 



The Gasser-Leutwyler counterterms also contribute to the cancellation of 
divergences in the expressions relevant to this paper. The subtraction is defined 
similarly to that of the weak counterterms, 



i 



16tt 2 



1 1/ 



1 - ln47r) 



r,, 



(21) 



with the divergent parts of the counterterm coefficients given in Table 2 |2Zj • 
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Table 2: The divergences in the strong 0(p 4 ) counterterms, r.-|27|. 



i 


r, 


1 


3/32 


2 


3/16 


3 





4 


1/8 


5 


3/8 


6 


11/144 


8 


5/48 



4 Partially Quenched Chiral Perturbation The- 
ory 

There are two approaches to (partially) quenched QCD, the supersymmetric 
formulation and the replica method [22]. Damgaard and Splittorff claim 
that the two methods are equivalent in the context of perturbation theory in 
the strong sector. We choose to follow the original method of Bernard and 
Golterman J5T for partially quenched chiral perturbation theory (PQChPT). 
In this method, the valence quarks are quenched by introducing "ghost" quarks 
which have the same mass and quantum numbers as the valence quarks but 
opposite statistics. As in |3()| . we consider a theory with n quarks and N 
sea quarks, so that there are n — N valence and n — N ghost quarks. The 
valence quarks have arbitrary mass, while the sea quarks are all degenerate. 
The symmetry group of the action is SU(n\n — N)l <8> SU(n\n — N)r. 
In the partially quenched case, the chiral field 



exp 



2i(j> a X a 



(22) 



/ 

has <fi a A a replaced by a (2n — N) x (2n — N) matrix, 

-(!$). <»> 

where 4> is an n x n matrix containing the pseudoscalar meson fields comprised 
of normal valence and sea quarks. (f> is an (n — N) x (n — N) matrix comprised of 
ghost-antighost quarks, while x is an n x (n — N) matrix of Goldstone fermions 
comprised of quarks and anti-ghosts. The most general set of operators can be 
constructed out of S, and these operators can be written in block form as 

where the sub-matrices have the same dimension as the elements of above. 
The transition to the partially quenched theory is made by replacing <f> a \ a by the 
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above (2n — N) x (2n — N) matrix, <I>, and replacing the traces in the operators 
with supertraces, defined as 

str(f/) = ti(A) - tr(D). (25) 

As a practical matter, in almost all of the NLO diagram calculations consid- 
ered in this paper, the minus sign in the supertrace is cancelled by an additional 
minus sign coming from anticommuting pseudo-fermion fields. The bare mass 
of a pseudoscalar meson is given by 

rriij = B Q (mi + rrij), (26) 

where rrii and rrij are the masses of the two quarks that form the meson. We 
define 17133 to be the tree- level meson mass of two valence strange quarks, as in 

123 

ml 3 = 2m 2 K -ml. (27) 

The tree-level mass of a meson made from the zth valence quark and a sea 
quark is 



m 2 s = S (to, + m s ) = -(ml + m| g ), 

i = u,d,s. (28) 

In this paper we consider only the case where the 7/ has been integrated out. 
Thus, the results are applicable to lattice calculations only when both sea and 
valence quark masses are small compared to m^i . Although it may be difficult 
computationally, this is precisely the case in which the LEC's of PQChPT are 
the same as those of full QCD when the number of sea quarks is three |20| . 
This is because the LEC's are independent of quark mass even if one varies sea 
and valence masses separately. In order that the LEC's of PQChPT be those 
of the real world, the sea and valence quarks must be small enough that the 
77' decouples, and its effects are integrated out the same way in both PQChPT 
and in full ChPT. 

The Minkowski space propagators for the flavor diagonal elements of $ are 
given by analytically continuing the Euclidean expression in |3(Jj 



13 p 2 - to^ + ie N \p 2 - m% + ie (p 2 - to| 4 + ie)(p 2 - m?- + ie) 



Ay = 
where 



Sijti 1 / 1 - m ss 



+1, for 1 < i < n (valence and sea); 
-1, for n + 1 < i < In- N (ghost). 



(29) 
(30) 
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At LO (NLO), the operators in PQChPT are still given by Eqs. (13,14) [Eqs. 
(16,20)], but with tr — + str for all operators. In the extension to the partially 
quenched case, 

Moo)' < 31 > 

in block diagonal form, and the mass matrix, 

X —> diag(m„, md,m s ,m sea , ...,m u ,m d ,m s ). (32) 

There is a choice in how to embed the quark charge matrix in the partially 
quenched theory, and this will affect the AI = 1/2, (8,8) amplitudes considered 
in this paper. If one wants to partially quench the electroweak penguin oper- 
ators, then the ghost quark charges should be the same as the corresponding 
valence quark charges. If, on the other hand, one wants to allow valence quarks 
to couple to photons and Z's, then the ghost quark charges should be set to 
zero so they do not appear in, and therefore cancel, the electroweak valence 
quark loops. We present amplitudes in this paper for both choices. Also, since 
we choose the sea quarks to have degenerate mass, the sum of the sea quark 
charges is the only quantity involving the sea quark charge that contributes. 
This is zero for three flavors, and in this paper we keep this true for arbitrary 
sea quark number, N, by setting the sea quark charge to zero. 

Also in the partially quenched case, the coefficient of the counterterm di- 
vergence depends on the number of sea quarks, N |23| . The N dependence of 
the necessary coefficients for the (8, l)'s was calculated following 28, 33], and 
the results are presented in Table 3. This paper uses a different basis from |2Sj 
for the (8, l)'s, and also several more LEC's appear here, so the calculation was 
redone for this work. The usual method was employed, expanding the action 
around the classical solution (background field method) and using a heat kernel 
expansion. The N dependence of the coefficients of the divergent parts of the 
(8,8) counterterms was given in ^Hj- These values are also presented in Table 
3. 

It is necessary to include an additional (8,1) operator, = str[Ae£ 2 ]str [5] , 

in this analysis of the partially quenched case since it can no longer be written as 
a linear combination of the other operators via the Cayley-Hamilton theorem. In 
the case of full ChPT, the operator is absorbed into the other operators, 

O 1 o' 1 ' ) , Cii' 1 ^, Ci2' 1 ' ) an d Ci3 1 ' ) - Since en has a divergent part, the coefficients 
of the divergences of the other four operators are modified (for N = 3) from the 
values in Table 1. Note that and ei2 have been omitted in Table 3. These 
LEC's do not appear in any of the amplitudes of interest in this paper. 

The Gasser-Leutwyler counterterms also contribute to the cancellation of 
divergences in this paper in the partially quenched case. The N dependence of 
the coefficients, Ti, is given in Table 4. 

When N is arbitrary, there is in general another operator tT\L il L u L^L v \, 
which cannot be absorbed into the first three Gasser-Leutwyler operators as 
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Table 3: The N dependence of the divergences in the NLO counterterms, e^'s 
and Ci's for the (8,l)'s and (8,8)'s, respectively. 
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Table 4: N dependence of the divergences in the strong 0(p 4 ) counterterms, 
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it can for N = 3 using trace relations. For the purposes of this paper, the 
additional operator and its divergent coefficient, Lq, can be absorbed into L\ 
through I/3 for the only amplitudes of interest to which it contributes, K — > nn 
for rriK = fn v , and in general, for UKX. Thus, we absorb the ./V dependence of 
Lq into Li,Li2 and L3 in Table 4. 

4.1 Role of the bilinear (3,3) operator 

The bilinear (3, 3) operator is useful in removing the power divergent coefficients 
to all orders in ChPT. Recall that the AI =1/2 matrix elements of the four- 
quark operators in general have a power divergent part. This power divergence 
reduces to a quark bilinear times a momentum independent coefficient 0- The 
quark bilinear operator can be defined as in 

Q^ = s(l- l5 )d (33) 

which is equal to a*- 3,3 ^Tr(A6S) to lowest order in chiral perturbation theory, 
where in our conventions, a( 3:3 ) = ^-B . As illustrated in Section 6, the matrix 
elements of this operator can be used to eliminate the power divergences in the 
effective four-quark operator matrix elements |34j . This subtraction is to all 
orders in ChPT, and in Section 6 we demonstrate this explicitly to NLO in the 
partially quenched theory, following the derivation in 0. It is crucial that the 
subtraction be independent of ChPT, since the higher order corrections of the 
power divergent operator can far exceed the physical contributions that one is 
trying to calculate. In order to carry out the argument to NLO in (PQ)ChPT for 
the case of the (8,l)'s we need the NLO LEC contribution of the O^ 3,3 ) operator. 
The effect of the subtraction involving this operator is to eliminate the LEC, 
CH2, at leading order, and to transform the NLO (8,1) coefficients to a linear 
combination involving the Gasser-Leutwyler coefficients. The chiral rotation 
eliminates the power divergent scale dependence (proportional to 02) of the 
LEC's to NLO. The effect of this transformation on the individual coefficients 
is given in Table 5 2 . 

5 K —> 7T7T without 3 momentum insertion 

In ^31 we have shown that all the amplitudes of interest for the (8,l)'s and 
(27,l)'s can be obtained to NLO in ChPT when one uses lattice computations 
from K° — > K°, K — > |0), K — > ir with momentum and K — ► tttt at the two 
unphysical kinematics points UK1 => mx = and UK2 JH| mK = 
Im-K- Specifically, these two points correspond to threshold, and, thereby, the 
Maiani- Testa theorem is evaded [0]. Here we ask how far one can get by not 
using 3-momentum insertion in K — > -k and using only non-degenerate quarks 
so that on the lattice m,K =/= m^. In this case one is using energy insertion 
with q 2 = (mpc — m^) 2 . The motivation for this should be clear. Not only can 

2 The table is constructed using information given in 1281 
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Table 5: When the tadpole terms are subtracted via the O^ 3,3 - 1 operator, the 
(8,1) NLO coefficients are transformed to new linear combinations involving the 
Gasser-Leutwyler coefficients. These new combinations no longer have power 
divergences. 

Transformed Coefficients 
el^e\-(Aa 2 /n(2Ll + H r 2 ) 

e r 2 ^e r 2 - (16a 2 /p)L r 6 
el^el + {Aa 2 /p)(-2Ll + Hl) 

e r io - elo ~ {^ 2 /f 2 )Ll 
eh - e r 13 - (8a 2 / f)L\ 



3-momentum insertion add to the computational cost, it also tends to be noisy. 
On the other hand, in a typical weak matrix element calculation, rriK ^ rn^ is 
relatively inexpensive to implement, since light quarks with several masses are 
needed anyway. 

At 0(p 4 ) in K — > tt one can see explicitly ^3] that different LEC's appear 
in front of (pk • P-k) 2 than in front of m^m 2 .. In general px • Pit 7^ mum-m so 
it is not clear if all of the LEC's needed for constructing K — > tttt to 0(p A ) can 
be obtained if one restricts to no 3-momentum insertion in K — > tt. We find 
that for all cases of interest without 3-momentum insertion, although some low- 
energy constants cannot be obtained, the linear combinations that are needed 
for constructing the physical NLO amplitude can always be obtained. This 
reduces the necessary computational effort considerably. This section will be 
restricted to demonstrating this result for the full theory, but in the next section 
we show that the same result holds also for the partially quenched case in the 
PQS framework. It is not known whether this continues to hold in the PQN 
framework, which is considerably more complicated at NLO. 

In ^21 we showed how to get physical K — ► tttt amplitudes for both AI =1/2 
and 3/2 cases to NLO. Since K — > tt amplitudes do not conserve four-momentum 
for m s rrid, it is necessary to allow the weak operator to transfer a four- 
momentum, q = px — P-m as in jl()| . This is also necessary for the case of 
K — > 7T7T at rriK = |14j . Our method |13j requires computation of K — » tttt 
at unphysical kinematics because there are low energy constants which appear 
in K — > tttt but do not appear in K — > tt at all [111 I3()j . 

There exist other unphysical kinematics values (besides UK1 and UK2) for 
the K — > tttt amplitudes where the initial and final state mesons are at rest 
that are accessible to lattice calculations. These values of the kinematics bypass 
the Maiani- Testa theorem because the final state pions are at threshold, but 
energy insertion (or removal) at the weak operator has to take place in order 
to conserve 4-momentum. These amplitudes have no imaginary parts as long 
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as nix > Tn n , (m s > m Uj 4, since then a two kaon intermediate state cannot go 
on-shell), and so bypass the Maiani- Testa theorem. On the lattice, so long as 
one studies the appropriate correlation function as a function of Euclidean times 
and does not sum over the time index, the weak operator can insert (or remove) 
the necessary amount of energy. What we are calling UK1 [tuk = m^) and UK2 
(mx = 2™^) are just special examples of this more general kinematics which 
we call UKX, which is itself a special case of the SPQcdR kinematics (one pion 
at rest, the other with 3- momentum inserted) [12j where both pion 3-momenta 
are zero, and E^, the energy of each pion, is equal to m,. 

We point out that the UKX kinematics is at threshold because of the ability 
of the weak operator to inject or remove the necessary energy, so that the 
Maiani- Testa theorem is bypassed even for AI =1/2 amplitudes. As pointed 
out by |15j . the case where wik — has a number of difficulties, especially in 
the partially quenched theory. It is, therefore, necessary to consider the more 
general kinematics of UKX (with mx > m^) in order to bypass this problem. 
Using UKX, one can then obtain all of the LEC's necessary to construct the 
(8,1), K — ¥ 7T7T amplitudes in both the full theory and in the partially quenched 
case, if a numerical calculation at UK1 is difficult or impossible. We present 
NLO results for UKX in the partially quenched theory in Section 8. 

Finally, it is also useful to emphasize that even when one works to LO, K it 
with nix ^ (without 3-momentum insertion) suffices to give K — ► nn at 
that order, thus providing an alternate subtraction method to the one that has 
been used recently [H El El d with K -► . 

5.1 (27, 1), AI = 3/2 

The expression for the physical K —> 7T7t, including only tree level 0(p 2 ) and 
0(p 4 ) weak counterterms, is ^2] 

(7T+7T-\O^ 7 ^^\K ) 



The counterterm expressions needed to construct this physical amplitude are 
given in ^3] Eqs (21-24), and the finite logarithmic contributions are given there 
in Appendix C. Counterterms needed to construct the above K — > nir amplitude 
can be obtained from K° — > K°; K + — * 7r + , AI = 3/2 (non-degenerate quarks); 
and K — > irir,AI = 3/2 at only one value of the unphysical kinematics (e.g. 
UK1 3 ). Note that the expression for K + — > 7r + , AI = 3/2 reduces, for the case 
of no 3-momentum insertion, i.e. q 2 — (wik ~ wi-^) 2 , to 

3 For the A7 = 3/2 case there is no difficulty at UK1. 
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+(^20 - dl~ 2d 7 )m 3 K m„ + (d 20 - d\ 
-d?)mKml], (35) 

The logarithmic corrections to this expression reduce to the value given in 
Appendix C of this paper. Fits to the K — > tt data can therefore give d 7 , 
d 2Q — d\ and d r 2 — 4^4 • Using these in the K — » tttt amplitude at the unphysical 
kinematics point (UK1) tuk = m n = m (Eq (23) of [13]) gives d\ — d r 5 . The four 
linear combinations [d 2 — ^d 24 , dj, d\ — d 7 5 , d\ — d 20 ] are sufficient to determine 
K — > tttt, AI = 3/2 at the physical kinematics as given in Eq.(34). Comparing 
Eq.(35) with the more general case of 3- momentum insertion, Eq.(22) of |13|. 
we see that the latter allows for separate determinations of d 2 and d 24 , whereas 
the simpler case of itik ^ without 3- momentum insertion, Eq.(35), gives 
only the linear combination d 2 — 4d 24 . Nevertheless, that suffices to get the job 
done. 



5.2 (8,1) + (27,1), AI = 1/2 

Recall that this is the most complicated case. The counterterms necessary to 
construct 0(p 4 ), [(8, l) + (27, 1)], AI = 1/2, K — > tttt amplitudes relevant for op- 

1/2 

erators such as Q 2 , which are mixed, can be obtained from the above values for 
e£'s and from the following AI = 1/2 processes: K° -> 0; K+ -> tt+ , AI = 1/2 
(non-degenerate quarks); and K — > tttt,AI — 1/2 at two unphysical kinemat- 
ics. All of the needed counterterm amplitudes appear in Section 4b of an d 
the corresponding logarithmic corrections appear in Appendix D of that paper. 
Note that an error was discovered since publication of that work in Eq (31) and 
in Appendix D, Eq (D6). The correct expressions appear here in Appendix F. 
Again, it is sufficient to allow q 2 = (rriK — m-ir) 2 in the expression for K — > tt, 
13 Eqs (28) and (29). These equations become 



+ (-6d[ - 2d r 2 + %d r 2i )m 2 K ml 

+ (-d r 20 + d\~ 3d r 6 + 2d r 7 )m 3 K m n 

+ (-d , 20 + d r 4 + 3d r 6 + d r 7 )m K ml], (36) 



+ {e r 2 + 2e r 3 + 2e r 5 - &e r 39 )m 2 K ml + {2e r 35 - 2e r 10 )m 3 K m n 
+(2e r 35 -e r n )m K ml\, (37) 

The logarithmic corrections associated with the above two amplitudes are 
given in Appendix C of this paper. In evaluating, for example, (ir + \Q 2 ^ 2 \K + ), 
the right hand sides of Eqs. (36) and (37) have to be added. In fitting to 
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lattice data, for example, the coefficient would give the combination [6d[ + 
2e\ + 2e r 2 - 2e£]. Also, in comparing Eq. (37) with Eq. (29) in [13] without 
3-momentum insertion, one can no longer separately obtain (e£ + 2eg + 2e^) and 
— 8e3 9 but only their sum; however, this is again sufficient to obtain the physical 
K -> irir amplitudes [Q3] Eqs (34), (35)] to NLO. 

We point out that for the AI = 1/2 amplitudes there are power divergences 
that must be subtracted using the 0( 3 ' 3 ) operator introduced at the end of 
Section 4. ft is crucial that the subtraction is to all orders in ChPT, since 
the higher order corrections of the power divergent operator can far exceed the 
physical contributions that one is trying to determine. This is discussed in more 
detail in Sections 7 and 8 for the partially quenched case, where we follow the 
derivation in given for the leading order case in the full theory (although 
there the analysis was done with quenched data) . The result of the subtraction 
is to eliminate the power divergent coefficient, a 2 , and to transform the (8,1) 
NLO LEC's to the values given in Table 5. Thus, fits to the subtracted lattice 
data will give the transformed coefficients, where their power divergences have 
been eliminated. This is what we want, since only these finite combinations 
appear in physical quantities. The process described in the above discussion on 
the determination of the NLO LEC's, along with that in J^jj is not invalidated. 

One can determine, using the AI = 3/2 amplitudes, the following constants: 
[dl,d 2 — 4^24,^7, d\ — dg,c?4 — d 20 }. H ere d\ an d dj can both be determined 
from K — > K [ 13 , Eq(21)], and the procedure for the others is given in the 
previous section. Given these, one can obtain e 2 rot and e\ rot — e£ rot from 
K° — > 0. Note that the values of the coefficients obtained are those of the 
subtracted amplitudes, and that the subscript refers to the LEC after the chiral 
rotation of Table 5 has been performed. Only after the subtraction can one fit 
to the lattice data using ChPT. Given the previous information one can obtain 

ei.rot + e 3,r°* ~ 4^9, An,rai ~ e 35 + l d 6> and 2e 10,rot ~ e ll + H from EqS 

(36) and (37), after the subtraction has been performed. From Eqs (30) and 
(31) 4 of reference ^3] f° r K — > nn, ttlk = rn v = m (UK1 5 ), one can then obtain 
e' ll +2e\ b rot — 'Sd r & . Making use of all of the input thus obtained into Eqs (32) and 
(33) of reference ^21 f° r K —> m K = 2ra. K (UK2), yields e\ 3 rot ~ |dg (after 
the subtraction). Thus, the 11 linear combinations necessary to construct the 
physical K — > 7T7t at NLO (without using 3-momentum insertion but with non- 
degenerate quarks in K — * w) are [d\ , d r 2 — 4e?2 4 , d^ , d\ — d\ , d\ — d r 2 Q , e 2 rot , e\ ro t + 
e 3.rot~ 4e3 9 , e-\Q tT0 t— 635+2^61 ^ e io,rot~ e ii+6^6' e n +2e^ 5 rot — 3c?g , e^ 3jrot — |dg]. 

5.3 (8,1) 

The case of pure (8,1) operators, e.g., Q 6 , is simpler than the previous case 
of mixed AI =1/2 operators, and is phenomenologically the most important 

4 Note that Eq (31) of 1131 is corrected in Appendix F, but this does not change the 
conclusion here. 

5 Although 1151 have pointed out that UK1 may be computationally demanding even for 
the full theory, it is not ruled out. In any case, for extracting the LEC's one can use the more 
general kinematics which we call UKX, as discussed earlier in this section. 
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one as it gives the dominant contribution to the CP-odd phase of e'/e com- 
ing from Q CD-penguins. For this case the six needed linear combinations are 

[ e 2,roti e l,rot + e 3,rot ~ ^ e 39' e 35 — e 10,rot' % e 35 — e\ x , e\ x + 2e r lbrot , e r 1Srot \. The 

first of these is obtained from K — > 0. The second requires both K — > and 
K — > tt (mx 7^ ^tt)- The third and fourth are also obtained from K — > 7T. 
X — > 7T7T at UK1 then gives the fifth, and K — > tttt at UK2 gives the sixth 
coefficient. Since it is likely that UK1 will prove to be particularly difficult |15| . 
it is possible to use another set of allowed values of UKX in order to obtain the 
remaining coefficients. Of course, one will want to do such a calculation using 
UKX anyway for the additional redundancy. All LEC's are those that would 
be obtained from a fit to lattice data after the power divergent subtraction has 
been performed. 

5.4 (8,8) 

Since the leading order (8,8) begins at O(p ), the NLO contribution comes at 
0(p 2 ). As an example, Eq (36) from [TO] is given (with our normalization of / 
and our convention for the Cj's), 



(^10(8,8)1^0 



2V2 



7C1 + c 2 



gC 3 I PK ' Pit — 



(38) 



Now with rriK 7^ fn^, even when both mesons are at rest, and pk • Pn = mK^m 
there is no loss of information, and all the coefficients can be obtained at NLO 
without 3-momentum insertion. 



6 Calculating K -> tttt Amplitudes in PQChPT 

In this section we discuss the ambiguity of PQChPT in the AI =1/2 case where 
eye-diagrams appear. At least two ways arise in the context of PQChPT for 
dealing with the gluonic penguins, the PQS and the PQN methods. These are 
described, and their predictions at leading order are compared using formulas 
given by Golterman and Pallante. In the following subsections we give NLO 
expressions in PQChPT for the ingredients necessary to obtain K — * tttt at 
0(p 2 ) and C»(p 4 ) for the (8,8)'s and (8,l)'s, respectively. For the (8,8)'s it is 
necessary to know K — ► tt, AI — 3/2 and 1/2 in order to get all the coefficients 
at NLO, as shown in [TO]. This remains true in PQChPT. The important point 
to note is that one can construct K — > tttt amplitudes for the (8,8) operator to 
NLO using only K — > tt with degenerate quark masses (uik — Tn^), along with 
K — > to perform the AI =1/2 power subtraction. 

For the (8,l)'s, one needs K — > 0, K — > tt with non-degenerate quarks, 
and K — > tttt at two values of unphysical kinematics, e.g. ttlk — m n (UK1) 
and rriK = 2m T: (UK2), as shown in ^3] in full ChPT to NLO. We have also 
introduced in Section 5 the kinematics for K — > tttt accessible to the lattice which 
we have called UKX, of which UK1 and UK2 are special cases. Reference ^H] 
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has demonstrated that UK1 has difficulties in the full theory, and is not tractable 
in the partially quenched theory due to enhanced finite volume effects. One can 
still obtain all of the needed LEC's to construct K — > tttt to NLO for the (8,l)'s 
from UKX, however. This remains true in PQChPT only when one is working 
within the PQS framework. This paper, therefore, follows the prescription of 
the PQS method for the (8,l)'s. Note that in the PQN method it is not clear 
if all the ingredients needed for constructing the physical K — > tttt amplitudes 
to NLO can be determined from the lattice, except in the full theory (N = 3, 
m sea — m va i) where the two methods coincide. Note, also, that the (8,1), 
K — > 7T7T amplitudes at UKX are afflicted by enhanced finite volume corrections 
except when m sea = m u = m<i for both the PQS and PQN methods. 

In the PQChPT case (as in full ChPT) the K — ► 7r amplitudes require non- 
degenerate quarks, m s ^ m u — rrid, in order to extract all of the necessary 
LEC's from them. Since this amplitude does not conserve four-momentum, for 
m s 7^ rrid the weak operator must transfer a four momentum q = pk — Pit ■ 
The conclusion of the previous section that three-momentum insertion is not 
essential holds also in the case of PQChPT. 

The diagrams to be evaluated for the NLO corrections are shown in Fig 
1. The topologies are unchanged from |13|. although additional pseudo-fcrmion 
ghost and sea meson fields propagate in the loops. The renormalization of the 
external legs via the strong interaction must be taken into account. 

6.1 The Treatment of Eye Graphs 

There is a subtlety concerning the AI —1/2 amplitudes in the partially quenched 
theory, and this has been discussed by Golterman and Pallante for the case of 
the gluonic penguins |21l 1301 l3l)] . What follows is a summary of their work. 
To illustrate the subtlety, we discuss the situation for the Qe gluonic penguin 
operator, given by 



The right part of this operator is a sum over light flavors, q = u,d, s, so in the 
full theory the right hand part is a flavor singlet under the symmetry group 
SU(3)fl. In the partially quenched theory one has at least two options. One 
may choose to sum over all the quarks, including sea and ghost in which case 
the right component of the operator transforms as a singlet under the extended 
symmetry group; therefore, this is called the PQS (partially quenched singlet) 
option. In the second option, one may choose to sum in Eq (39) over only the 
valence quarks. In this case the operator is a linear combination of two terms, 
one of which transforms as a singlet under the extended symmetry group, while 
the other does not transform as a singlet under the irreducible representation 
of the extended symmetry group (rather, for Qg, it transforms in the adjoint 
representation); therefore, we choose to call this the PQN (partially quenched 
non-singlet) method. 




(39) 
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Figure 1: Diagrams needed to evaluate the NLO amplitudes in (PQ)ChPT. 
NLO corrections include tree-level diagrams with insertion of the NLO weak 
vertices (crossed circles), tree- level diagrams with insertion of 0(p 4 ) strong ver- 
tices (lightly shaded circles), one- loop diagrams with insertions of the LO weak 
vertices (small filled circles) and the 0(p 2 ) strong vertices (big filled circles). 
The lines represent the propagators of mesons comprised of valence, ghost, and 
sea quarks. Al and A2 are for K — ► 0. B1-B3 arc for K — > ir. C1-C6 and 
D1-D6 are for K — ► irir. 
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Given that the flavor blind, vector character of the quark-quark-gluon ele- 
mentary interaction in QCD plays a crucial role in leading to the explicit singlet 
form [Eq (39)] of the right-hand part of the penguin operator, it seems reason- 
able to preserve this basic character in generalization to the partially quenched 
case which contains additional quarks. This provides the rationale for the PQS 
option. 

The origin of the PQN option is quite different; it is, in fact, the straight- 
forward implementation of the quenched approximation to a lattice calculation 
of the necessary Green's functions. The usual practice leads one to use only 
the valence quarks in the necessary Wick contractions for, say (tt\Qq\K), which 
then lead to valence quark loops (see Fig 2a, b), the so-called eye graphs. In such 
an implementation all other quark loops are computed when the fermion deter- 
minant is evaluated in the generation of the gauge configurations. When one 
partially quenches in the PQN method, the gauge configurations are generated 
using the number and mass of the sea quarks, but the propagators for the loops 
of the eye graphs (Fig 2a, b) are still computed with those of the valence quarks. 
In the partially quenched case where the sum in Eq (39) is over the valence 
quarks only, as mentioned above, the operator is a linear combination of two 
terms, only one of which transforms as a singlet under the extended symmetry 
group. 

Fig 2 shows the Green's function relevant for a lattice evaluation of (tt\Qq\K) 
consisting of the two eye graphs originating from the Wick contractions. Any 
number of gluon lines from the background gauge configurations (not explic- 
itly shown) are understood in such a pictorial representation of these non- 
perturbative graphs. As usual, one of the Wick contractions is a product of 
two color traces (Fig 2a), while the second is a single trace over color indices. 
In the PQS implementation of the Qq penguin operator, in the quenched case 
where qq loops in the gluon propagation are not included, the eye graph (Fig 
2b) with a single color trace should also be excluded, for consistency [21 EH] ■ 

In the PQN option of calculating (tt\Qq\K), one uses valence quarks for the 
propagators of the eye-graphs in the corresponding Green's function, as this 
appears analogous to the usual practice in lattice computations. However, the 
situation at hand demands caution. Lattice calculation of (tt\Qq\K) is quali- 
tatively different in important aspects from (say) spectrum, decay-constant or 
form-factor calculations. 

To trace the potential inconsistency we show the weak operator with a mag- 
nified view in the non-perturbative eye-graph (Fig 3). Inside the dashed lines is 
the magnified short distance effective penguin operator; outside of these dashed 
lines any number of soft gluon lines from the background gauge configurations 
are understood, just as in Fig 2. For (it\Qq\K), Fig 3a and 3b correspond to 
the product of two color traces (Fig 2a) and Fig 3c corresponds to the single 
trace over color indices (Fig 2b). Fig 3c shows clearly that the corresponding 
Wick contraction (single trace over color indices for Qq, i.e., Fig 2b) in a lat- 
tice evaluation of (tt\Qq\K) contains a qq loop in the propagation of the gluon, 
and since in the quenched case these are being dropped from the background 
gauge configurations one may wish to exclude Fig 2b (for Qq) in the quenched 
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(a) 



(b) 



Figure 2: The quark contractions needed for K — » it, AI =1/2 matrix elements 
include the above eye diagrams. A connected line represents a trace over color 
indices, so Fig 2a represents a product of two color traces, whereas Fig 2b 
represents a single color trace. 

approximation. In a similar vein, for the partially quenched case one may, for 
consistency, take the quark loop in the eye graph of Fig 2b (Fig 3c) to be that of 
sea quarks only j^EH], as in the PQS method. This lack of consistent (partial) 
quenching causes the low energy dynamics of (P)QChPT to change between 
PQS and PQN methods. It is not clear if the additional (partially) quenched 
non-singlet terms that modify the low energy dynamics correctly account for the 
otherwise neglected loop contractions, or if the (partially) quenched low energy 
constants from the singlet operator alone substituted into the full ChPT formu- 
las for K — > 7T7T provide a better estimate for the physical amplitudes. Thus, 
the appearance of eye-diagrams has created an ambiguity because the contrac- 
tion of Fig 3c yields a quark vacuum bubble, and it is not obvious whether the 
propagators to be contracted should be the sea or the valence; again, the first 
choice corresponds to PQS and the second to PQN. 

The correspondence between the traditional form of the non-perturbative eye 
graphs as shown in Fig 2 and the non-perturbative eye graphs with the magnified 
view of the penguin operator as shown in Fig 3 for all penguin operators is as 
follows. For the Q3 and Q5 operators, the color contraction of Fig 2a corresponds 
to Fig 3c, while the color contraction of Fig 2b corresponds to Figs 3a and 3b. 
For the Q4 and Qq operators, the color contraction of Fig 2a corresponds to 
Figs 3a and 3b, while the color contraction of Fig 2b corresponds to Fig 3c. For 
the electroweak penguins, the picture in Fig 3 carries over, but with the gluons 
replaced by a photon or a Z. In that case, for Q7 the color contraction of Fig 
2a corresponds to Fig 3c, while the color contraction of Fig 2b corresponds to 
Figs 3a and 3b. For Qg, the color contraction of Fig 2a corresponds to Figs 3a 
and 3b, while the color contraction of Fig 2b corresponds to Fig 3c. In short, 
Fig 2a corresponds to Fig 3c for the operator Qi, i = 3 — 8, i odd, while Fig 2b 
corresponds to Fig 3c for i even. 

The treatment of ChPT for the case when only valence quarks are contracted 
in the eye-diagrams was first discussed by , for the case of the gluonic pen- 
guins. When one includes only the valence propagators in the eye-diagrams (no 
partial quenching of the effective operator) for the case of the gluonic penguins 



23 



u 

u 

(a) /' ( b ) 




LI 



(c) 

Figure 3: The quark contractions needed for K — ► 7r, AI = 1/2 matrix elements 
include the above eye diagrams. The weak operator is shown with a magnified 
view inside the dashed lines so that one can see how it arises in perturbation 
theory for the gluonic penguins. Fig 3a corresponds to q being contracted with 
d, while Fig 3b corresponds to q being contracted with s. Fig 3c corresponds to 
q being contracted with q. For the electroweak penguins, one would replace the 
gluon lines with those of photons or Z's. 

the right hand part of the (8,l)'s is no longer a singlet, and there is a contribu- 
tion from a non-singlet operator. For the left-left gluonic penguins, Q3 and Q4, 
these non-singlet contributions do not occur until next-to-leading order |36j . For 
the left-right gluonic penguins, Q5 and Qq, the non-singlet operator transforms 
under the same irreducible representation as the (8,8) electroweak penguins. 
Since the (8,8)'s are NLO at 0(p 2 ), even the leading order gluonic penguins 
can have logarithmic contributions from the one loop insertions of the lowest 
order (8,8) operator. These were calculated in [50J for the left-right gluonic pen- 
guins, Q5 and Qq. Since the amplitudes in this case no longer transform as pure 
(8,l)'s, but pick up a contribution from the (8,8)'s, this calculation corresponds 
to the PQN method. It is useful to compare the (PQ)ChPT expressions for the 
PQS and PQN methods, and the next section compares the two methods at 
leading order for the left-right gluonic penguins, using expressions derived by 
Golterman and Pallante [2T1 1501 150| . 
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We choose to work within the framework of 130], where the PQS method 
was (implicitly) used. In this case there is the possibility of determining the 
LEC's to NLO. For the left-right gluonic penguins, for example, at NLO in the 
PQN method there are many more LEC's that appear in the amplitudes we are 
considering than in the PQS method. These are the 0(p 4 ) LEC's of the (8,8) 
NNLO local operators, and it is not even clear whether one can determine the 
correct linear combinations of the new LEC's necessary to construct K —* tttt 
at NLO in ChPT from the PQN method, except when N = 3 and m sea = m va i 
(i.e., full QCD), as in that special situation the two options coincide. On the 
other hand, for the PQS method, no new ingredients are needed over the ones 
listed in our previous work J3j which were needed for the case of full ChPT 
6 . The PQS method can also be applied to obtain all of the needed LEC's to 
construct K — > tttt to NLO for the case N = 2 (using the same ingredients as 
for the full theory), though in this case the LEC's are not necessarily the same 
as in the N = 3 physical case. 

To reiterate, in general, the PQN method is complicated by the contributions 
of many more LEC's, and it is not known whether this method can be used to 
NLO. Such a determination would require a two-loop calculation. The PQS 
method gives us everything we need, and is the only method where we have 
demonstrated that it is possible to obtain K — > tttt to NLO in ChPT. Thus, we 
use the PQS prescription. 

As discussed in |5Tj, the NLO (8,8) LEC's that appear in linear combinations 
with the LO (8,1) LEC's in the PQN expressions for the amplitudes of the left- 
right gluonic penguins are not present for the PQS method. For the case where 
N = 3, there is no ambiguity, and one must extract the (8,1) LEC's separately 
since these LEC's take the same values as in the full theory. However, when N 
is not equal to 3, it may be that the additional (8,8) LEC's appearing in linear 
combinations with the (8,1) LEC's bring the N ^ 3 values of the (8,1) LEC's 
to closer agreement with the N = 3 values of the real world. As long as an 
explicit N = 3 lattice calculation is lacking, it may be useful to compare the 
determinations of both PQN and PQS leading order LEC's at other values of 
N in order to learn something of the size of the systematic error due to partial 
quenching [^J|31)|. This is discussed further for Qs,6, LO K — ► tttt amplitudes 
in the next subsection. 

This paper requires K — > 0, K — ► tt, m s ^ m<j = m u , and K — > tttt at two 
unphysical kinematics in order to construct the physical (8,1), K — ► tttt am- 
plitude for any gluonic penguin operator (^3,4,5,5) using the PQS prescription. 
Reference j3U] presented K — > 0, and K — » tt, m s — m u — ma, and we agree 
with those calculations in the case we consider, namely the partially quenched 
case with m va i , m sea -C m„> . We extend these calculations to include all am- 
plitudes needed to obtain the (8,1) LEC's necessary to construct K — > tttt to 

6 Note that 1151 have shown that there are difficulties at what we call UK1 (K — > tttt 
with m,K = m^). There is, however, an additional set of kinematics points that bypass the 
Maiani- Testa theorem, creating the two pion state at threshold with energy carried by the 
weak operator which we call UKX. If UK1 proves difficult or impossible one must supplement 
the other ingredients with a calculation of K — > tttt at UKX with m,K > m n . 
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NLO. 

In the case of the (8,8), AI =1/2 amplitudes, one must also make this 
choice of whether to (partially) quench the right side of the penguin operator. 
In this case, however, the difference comes in the choice of the quark charge 
matrix, Q. If we choose the ghost quark charges to be equal to the valence 
quark charges, then we quench the electroweak penguins, and one should ignore 
the valence contributions to Fig 3c (with the gluon replaced by a photon or Z, 
Fig 3c corresponds to Fig 2a for Q 7 and to Fig 2b for Qg) in the lattice calcu- 
lation. However, if one chooses the ghost quarks to have zero charge then the 
electroweak interaction remains unquenched, and one must include the valence 
quarks in Fig 3c in the lattice calculation. In both cases the sea quark loop con- 
tributions to Fig 3c vanish if we assume degenerate sea quark masses and that 
the sum of the sea quark charges is zero. The logarithmic expressions resulting 
from either choice for the (8,8) 's are presented in Appendix D. 

To summarize, our calculation for the (8,1) gluonic penguin matrix elements 
corresponds to the PQS method. In the corresponding lattice calculation, the 
eye contractions of Fig 3c (corresponding to Fig 2a for Q 3 , Q 5 , and Qr, and to 
Fig 2b for Q4, Q 6 and Q$) include only the sea quarks. That is, the propagator of 
the internal loop of Fig 3c is calculated with the masses of the sea quarks, not the 
valence quarks. As discussed above, when N = 3, this, the PQS method, allows 
for the only known implementation of the reduction method for the gluonic 
penguins. It greatly simplifies the LO analysis and makes possible a NLO 
determination of all of the necessary LEC's, as demonstrated in this paper. For 
the (8,8) electroweak penguin matrix elements, for degenerate sea quark masses 
the eye graph of Fig 3c (with the gluon replaced by a photon) vanishes for any 
number of dynamical flavors by construction (see Sect 4). Whether one chooses 
to include valence quarks in the loop of Fig 3c does not significantly alter the 
situation in PQChPT, and formulas for both implementations are given in this 
paper. 

6.2 PQS vs PQN at Leading Order 

This section is a review of Golterman and Pallante's |211 131)1 13f)| results for the 
leading order, left-right gluonic penguins, Q5 and Qq. Table 6 compares the 
results of the PQS method versus those of the PQN method. The results are 
for the subtracted K — > 7r matrix elements, where the (large) subtraction is 
performed using K — > 0. For details on how this subtraction is performed, see 
0. The end result of this subtraction in the case of full QCD (no quenching) 
is just a^ =3 , which is the physical LO LEC that contributes to K — > hit. In 
this case, the two methods, PQS and PQN, are procedurally the same, and they 
therefore give the same answer. 

When N = 3, but m sea ^ m va i, the LEC's in the amplitudes are still those 
of the full theory, but an additional LEC, the leading order (8,8) electroweak 
penguin LEC, ass, contributes in the PQN case to K — > multiplied by some 
logarithmic terms [23- Thus, a subtraction that is performed without taking 
this into account has a contamination. That is, there is an extra term appearing 
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Table 6: The leading order LEC's as determined in PQChPT from K — ► ir after 
using the K — > subtraction described in are presented. They are compared 
for the PQS and PQN methods in the case of the left-right gluonic penguins. 
The two methods agree for the full QCD case. For the N = 3, m sea ^ m va i 
case there is a logarithmic contamination for the case of PQN. That is, there 
is an extra term appearing at leading order that must be accounted for in fits 
used to obtain the subtraction coefficient from K — > 0. For N = 2, the LEC's 
are not those of the full theory, and additional terms appear for the PQN case. 
For N — 0, the quenched case, there are also additional terms that contribute 
in the PQN case. See j^OE] for the derivations of these results and the values 
of the logarithmic corrections abbreviated here. 





N = 3, m sea = m val (Full QCD) 


PQS 


PQN 


<= 3 






N = 3, m sea ^ m va i 


PQS 


PQN 


<= 3 


a N=3 _|_ Q,N=3(/ g terms) 




N = 2 


PQS 


PQN 


<= 2 


K =2 + W^(A S ' S) + |/f ' 8) ) + a^2(/o 5 terms) 




N = 


PQS 


PQN 




K = ° I^F^f S + 5^ S ) + «Q 5 (^5 terms) 
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at leading order that must be accounted for in fits used to obtain the subtraction 
coefficient from K — ► 0. Looking in Table 6 at the PQS result, we see that this 
method is simpler. At NLO the difference is even more severe, so that PQS is 
the only method shown to be feasible. In this case the difference in practice 
between the two methods is whether one uses the sea mass or the valence mass 
in the propagator of the loop in Fig 3c. (See the preceding section for the 
correspondence between Fig 3c and the traditional form of the eye-diagrams in 
Fig 2 for the various operators.) 

When N = 2 the LEC's are no longer those of the full theory, and an 
ambiguity results. In this case the calculations differ in that for the loop of 
Fig 3c for PQS one uses 2 flavors with the dynamical mass, while for PQN 
one uses 3 flavors with the valence masses. Here, a|^ =2 appears multiplied by 
logarithmic terms, and these must be removed in the fits to K — > before 
the subtraction can be performed. Notice also the presence of the /3 terms in 
linear combination with the a^~ 2 term. These j3 terms always appear in the 
same linear combination with a^ =2 , including in the expression for K — > itit. 
Thus, it is not obvious whether they represent a correction to the a^~ 2 term 
or a contamination. Clearly, it will be important to compare the results of 
both methods. Note also that the f3 terms that appear with af =2 have a scale 
dependence proportional to a^ =2 , and that if a^ =2 is not so far from a^ =0 , as 
determined in [2], then this scale dependence would be large. It would then be 
necessary to include the partially quenched chiral logs proportional to a% s =2 in 
K — ► 7T7T in order to cancel the scale dependence and obtain a consistent answer. 
It would be extremely useful to have (at least so long as an N = 3 calculation 
is not available) a study of the subtracted LO constants for the PQS and PQN 
methods as a function of N, so that one could try to extrapolate each result to 
N = 3, and compare the two. 

Note that the (3 terms are related to the (8,8) LEC's, Cj in the terminology 
of this paper. The correspondence to our notation is 

M 8 < 8) = (4 7 r) 2 2 C3 , 
4 8 ' 8) = (4tt) 2 2 Ci , 

^ 8 ' 8) = (47r) 2 2c 4 . (40) 

Finally, when N = the theory is completely quenched. This corresponds to 
ignoring all contractions of the kind in Fig 3c for the PQS method and keeping 
them with the valence quarks for PQN. In the quenched case, the additional 
LEC's for PQN, the cxq S and (3f s , are coefficients of non-singlet operators, but 
they have no relation to the (8,8) electroweak penguins. Golterman and Pallante 
|2*T] provide a possible recipe for determining ag S on the lattice. Analogous 
to the partially quenched case, if ctg S is large, it would imply a large scale 
dependence on the subtracted combination of LEC's in the PQN method, which 
would have to be cancelled by the quenched logs proportional to a% s in K — > 7T7T 
in order to obtain a consistent, scale independent answer. There are indications 
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from the large N c (N c is the number of colors) approximation that the LEC, 
ocq S , is indeed large compared to =0 |2Z|. Again, it is not clear if the non- 
singlet terms represent a correction or a contamination, and results for both 
methods should be compared as part of an extrapolation in N . 



7 (Partially) Quenched (8,8) 's to NLO 

In this section we present the results for the partially quenched K ^ tt and 
K — ► amplitudes needed to construct the K — > tttt amplitudes to NLO for the 
(8,8)'s. The power divergent subtraction is discussed for the AI = 1/2, K — ► tt 
amplitude. Formulas are presented for K — » and K — > w for nondegenerate 
quark masses, as well as K — > tt for degenerate quark masses. It is demonstrated 
that K — > tt with degenerate masses is sufficient to construct K — > 7T7t to NLO in 
the partially quenched theory, while K —> n with non-degenerate quark masses 
gives additional redundancy in determining the NLO LEC's. 

We show that in the case of K — ► 7r with degenerate quark mass, the iV = 
limit of our expressions produces the quenched result, which will be useful for 
fits to already existing lattice data. It is important to notice that not all LEC's 
needed for NLO K — ► tttt can be determined from the quenched K — ► tt data, 
since Cg, which is needed in the physical K inr expressions, does not appear 
in the quenched K — > tt formulas. One can see the scale dependence of the 
LEC's from the formula, 

Ci{H2) = Ciim) + In—, (41) 

(4tt/) 2 ^ 2 

which can be obtained from Eq (19), the definition of the renormalized LEC's. 
The coefficients, r/i, are given in Table 3. Since the scale dependence of the Cg 
coefficient in the physical K — > tttt amplitude is significant (where it is needed 
to cancel the corresponding scale dependence in the NLO log terms) , it is crucial 
to do dynamical simulations of the (8,8) K — > n amplitudes in order to bring 
under control the systematic errors due to the chiral expansion. 



7.1 Partially Quenched (8,8)'s with nondegenerate quark 
masses 

The LEC's needed to construct the K -> tttt, AI = 1/2 and 3/2 (8,8) 's can be 
obtained from the K —t tt amplitudes with energy insertion and m s ^ rrid = m u . 
The (8,8) K — ► tttt counterterm contributions for both the AI — 3/2 and 1/2 
amplitudes are given by 



(7r+ ^ |O (8,8),(3/ 2)|x )ct = _^8| + « [{ - C r_ c r_ 2c r_ 2cl _ 4c r )m 2 

-{-c\-^+Acl + Acl + 2cl)ml], (42) 



29 



(*+*-|O<M>.(VS>|X0) et = -^--^[(-^-eS + ^ + ^ + ScSJm^ 

+ (-<£ + ^ + 2cJ + 8< + 8c£ + 44)mJ]. (43) 

These are the expressions in the full theory and were given by where 
they showed that one can obtain the necessary linear combinations of LEC's 
from K — > 7r with momentum, AI — 1/2, 3/2. We demonstrate this holds also 
for the partially quenched case (without the need for 3-momentum insertion, as 
explained in Section 5). In Eqs (42), (43) as well as all the following amplitudes, 
we include only the tree level weak counterterm contributions. For clarity, the 
logarithmic terms and the Gasser-Leutwyler Li counterterms have been omitted 
from this section, but are included in Appendix D. 

The K — > 7r counterterm amplitudes are given by 

(7r+ |0(8 > 8),(3/2) |x+)ct = ^ + j_ [2{cl + c r )m 2 K+2{cl + c r )m 2 

-{c[ + c r 2 )m K m„ + 2c r 6 Nm 2 ss ], (44) 



- (ci - c 7 2 - 2(Qm K m n + Ac 7 6 Nm 2 ss ] . (45) 

At this point, a practical issue in the extraction of the LEC's should be 
mentioned. There is a power divergence in the NLO coefficient c\ due to mixing 
with unphysical lower dimensional operators that must be removed if one is to 
have any hope of numerically extracting any of the LEC's. This is a problem for 
K — > 7r, AI = 1/2, but not K — > n, AI = 3/2, since the combination c\ + c r b is 
finite in the continuum limit. For the AI =1/2 amplitude the power subtraction 
method of RBC 2 can be used, and this requires the K — > amplitude. At 
NLO, this amplitude is 



(010(8,8)^0) = ^ [2Mm 2 K) _ Mml) _ Mm 2 3) 

+NA (m 2 sS ) - NA (m 2 uS )} - jcl(m 2 K - mj), (46) 

where Ao(m 2 ) is defined in Appendix A. We also mention that the (8,8), K — > 
calculation has an eye-diagram, and one must make a decision whether to keep 
the valence quarks in the eye contractions or not. The above formula, Eq (46), 
corresponds to keeping the valence quarks in the eye contraction. If one neglects 
the type of contraction associated with Fig 3c, then one obtains 
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(O|0W9|*O) = ^[NA ( m l s )-NA a ( m l s )]- S jcl{ m \- m l). 

(47) 

Unlike the case of Qq, the chiral perturbation theory is not substantially changed, 
and one can use either method, as long as one is consistent. The K ^> sub- 
traction works as follows. We make use of the subtraction operator introduced 
in Section 4, 

9 (3 < 3) = 3(1 - 75 )rf = a (3 ' 3) Tr(A 6 £) (48) 

to lowest order in chiral perturbation theory. The mass dependence of the above 
quark bilinear operator, 0( 3 .3) , is the same as that of the power divergent part 
of the four-quark operators, so one can use the matrix elements of this bilinear 
operator to subtract out power divergences to all orders in ChPT. In order to 
perform the subtraction at NLO for the (8,8) 's we need the following leading 
order expressions of the 6( 3 > 3 ) amplitudes, 

(7r +| (3,3)| ir+) = _2 a (3,3) j (4Q) 

(0\e^\K°) = ja^l (50) 
When we take the ratio of (0\O^\K°) to (O0 3 ^\K°) we get 

v ' _ 1 '- =-4— i =-(mi-r<) + 2— ^=-(logs) + ... (51) 

where we have omitted terms of higher order in the chiral expansion. Note, 
however, that all higher order terms proportional to c\ cancel in the ratio. 
Fitting to this expression allows one to obtain c r Ja^^\ which one can then use 
in the subtraction of the power divergences of K — > it. Notice that c\ has a 
scale dependence that must cancel the scale dependence of the ass log term in 
K — > 0. Thus, one must ensure the value of /i in a chiral fit to K — > n is the 
same as the value of \x used in the power subtraction. After the subtraction, 
the following expression no longer has power divergences. 



(7r+|o( 8,8),(l/2)| X+) + 4 ^ (7r+ |e(3,3)| X+) = ^88 (1 + /ogs) + ^ [4K + cD 

+4« + c r 5 )ml - (c[ -c r 2 - 2c r 3 )m K m v + ANc r 6 m 2 ss ]. 

Here, c\ appears only in the linear combination c\ + eg, which does not contain 
power divergences. Thus, the K — > subtraction has removed the power diver- 
gences from the AI = 1/2, K ^ it expression, including all of the higher order 
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power divergent contributions, an important point, since the subtraction does 
not require (PQ)ChPT for its implementation. The expression, Eq (52), is the 
one which should be fitted for the NLO LEC's. Thus, fitting to (44) and the 
power subtracted amplitude, (52), one can obtain all of the linear combinations 
needed for K — > irir at NLO. 

In principle, one can obtain a 88 from either leading order term. In practice, 
it is safer to use the 3/2 amplitude since the 1/2 amplitude could receive some 
residual chiral symmetry breaking contribution unless one uses a discretization 
that has exact chiral symmetry. The AI = 3/2 expression does not involve 
power divergent subtractions, and is, therefore, the best way to get the lead- 
ing order coefficient. One can get Cg from the term that depends on the sea 
meson mass. From fits to the other mass combinations, one obtains c\ + c 5 ', 
c\ + C2, and c\ — c!^ — 2cg. Along with ass, the four linear combinations: 
[c[ + C2,c\ — — 203,04 + Cg,Cg] are sufficient to determine K — > nn at the 
physical kinematics, as one can verify with some simple algebra from Eqs (42) 
and (43). When TV = 3, the values of the LEC's determined from PQChPT 
are the same as in the full theory. We point out in the next subsection that 
one can get all of the needed information to construct the EWP matrix element 
for K — ► 7T7T to NLO even with K — > 7r using degenerate valence quark masses, 
along with K — > to perform the power subtraction in the AI — 1/2 case. 
The nondegenerate case remains useful, however, in that it provides additional 
redundancy in determining the NLO LEC's. 

Note that for the cases of physical K — > mr (8,8) amplitudes (42), (43), and 
(58) for the corresponding (8,l)'s, the pseudoscalar decay constants and masses 
are the physical (renormalized to one-loop order) ones. For all other amplitudes 
given in this paper except K — > irir at physical kinematics, the formulas are 
in terms of the bare constants. The distinction between bare and renormalized 
constants is made only in tree-level amplitudes, since making this distinction 
in the NLO expressions introduces corrections at higher order (NNLO) than is 
considered here. 

The logarithmic and Gasser-Leutwyler counterterm contributions to the am- 
plitudes in this section are given in Appendix D. 

7.2 Partially Quenched (8,8) 's with degenerate quark masses 
for K — > 7r 

For the case of degenerate quark masses, Eqs (44) and (52) become 



( 7r +jO(8.8).(3/2)| A -+) 



4as8 



+ -^[(-cj - c r 2 + Ac\ + 4cg)m 2 + 2cl Nm 2 ss ], 

(53) 



ct 



p 




\K + ) ct 



8a 8 8 



+ jzK-cl + c r 2 + 2cl + 8c r 4 + 8c r 5 )m 2 + 4c r 6 Nm 2 ss }. 

(54) 



P 
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where the logarithmic parts of the above expressions are given in Appendix D. 
The K — > subtraction is performed exactly as in the non-degenerate case, 
yielding the above result, Eq (54). Again, simple algebra will verify that the 
above linear combinations of LEC's are sufficient to determine the LEC combi- 
nations in Eqs (42) and (43) for the physical (8,8) K — > irn amplitudes to NLO. 
For example, if one subtracts the m 2 coefficient in Eq (54) from the to 2 coef- 
ficient in Eq (53), one gets the same linear combination as the first four terms 
in Eq (42). If one wants to obtain all the needed information in the full theory, 
one must do an N — 3 simulation, varying the sea quark mass with respect to 
the valence quark mass in order to determine Cg. That is, one must still vary 
the sea quark mass independently of the valence quark mass. 

7.3 Quenched (8,8) 's with degenerate quark masses for 

K -> 7T 

One can obtain results in the quenched theory by taking the N — limit of the 
K — > and degenerate K — > tt formulas. One then obtains, for K — > 0, 

(O|0W»|ff°) = ^l[2A (ml) - A Q (ml) - A (mt 3 )] - jcl(m 2 K - ml). 

(55) 

In the quenched theory, the scale dependence of c\ vanishes, as one can verify 
from Table 3. Therefore, the scale dependence of the logarithms in quenched 
K — ► must also vanish; the fact that it does can be seen from Eq (55). The 
subtraction is performed the same way as in the partially quenched theory, and 
the expressions for K — > 7r are 

( 7r +| (8,8),(3/2)| A -+) 



4^88 

~F . 

4m 2 
+ — 



1 



167T 2 / 2 

-16a 8 8 , 



2 , m 2 
m In —5- + m 



cl-c r 2 + Ac\ + 44 



(56) 
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Vt -c r 1 + c r 2 



2cl 



8c\ 



-8c r 5 



(57) 



Here, is the quenched Gasser-Leutwyler coefficient that appears in f v . No- 
tice that Cg does not appear in the quenched theory, though it does appear in 



33 



K — > 7T7T, where, as one can see from Table 3 it has a non- vanishing scale de- 
pendence. This dependence on the chiral scale leads to a large uncertainty in 
the NLO K — > tttt amplitudes, since the scale dependence of the LEC's must 
cancel against those of the chiral logarithms. Thus, it is essential to compute 
K — > 7r with dynamical quarks in order to reduce the uncertainty due to the 
chiral expansion. 



8 Partially Quenched (8,l)'s to NLO 

The amplitudes necessary to construct the physical K — > tttt matrix elements 
are K — > 0; K — > n, m s ^ to^ = m u ; and K — > tttt at the unphysical kinematics 
points of UKX, of which UK1 and UK2 are special cases. The counterterm part 
of the physical (8,1) amplitude is 



(ejo - 2e^ 3 + 2e r u + e r 15 )m 2 K + (-2eJ + 2e r w 

+e r n + Ae r 13 + e r u - 4e r 2 - 2e\ - Ae r 35 + 8e r 39 )ml 
+ ^[2m 2 K L 4 + (-4L 4 - L 5 + 8L 6 + AL 8 )ml] 

(58) 

This differs from our previous expression [T3|, Eq. (35) in the appearance 
of a new LEC, e\ i: and in the inclusion of the Gasser-Leutwyler coefficients of 
the amplitude previously given separately as part of the log terms in D10 of 
|13| . The Gasser-Leutwyler coefficients are included here with the rest of the 
LEC's for clarity. In the full theory the operator corresponding to this LEC 
can be absorbed into the other operators of^, 0^{ X \ O 1 ^ and O^ 1 ' via 
the Cayley-Hamilton theorem, as discussed in Section 4. Since this is no longer 
true in the partially quenched theory, one must obtain the constant separately. 
Therefore, it is left explicit in the physical amplitude. 

For K° — > 0, we have 



<0|O (8 ' 1) |^°) 



4ia 2 , 2 2n 
— — [m K -mj 



-e r 2 Nm 2 ss ] 



j(m 2 K - m l)[(2e[-2el)m 2 K 



(59) 



The expression for K — > 7r is 



—airriKm^ 



2 Tc\/ r r\ 4 

j^a 2 m K - j 2 -[2{e 1 - e 5 )m K 
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-2(eJ - e r 35 )m 3 K m n + (2e r 3 + 2e\ - 8e r 39 )m 2 K ml 
+(2eg 5 - e r n )m K ml + Ne r 2 m 2 K m 2 ss - Ne\ A mKm,T,m 2 ss \. 

(60) 

The (8,1) amplitudes have power divergent parts that must be subtracted, 
just as in the case of the AI = 1/2, (8,8) 's. The subtraction is performed in the 
same way, but in this case the power divergent coefficient, Q2, is present already 
at leading order. Thus, we must consider the effects of the power subtraction at 
NLO if we are interested in obtaining the LEC's to this order. As we will see, 
the effect is to modify the NLO, (8,1) LEC's by adding terms proportional to 
the Gasser-Leutwyler coefficients. The LEC's modified in this way are just those 
that have a scale dependent part proportional to a 2 , and the new constants so 
obtained are given in Table 5. The power subtraction eliminates the tadpole 
contributions from the amplitudes, and the new combinations of LEC's in Table 
5 are free of power divergences, and can be obtained in numerical fits to lattice 
data. 

The ratio of the K -» amplitude to the 6 (3 ' 3) K -» amplitude to NLO 
in PQChPT is 



(0|Q(M)|jrQ) 

(010(3.3)1^0) 



Qf2 



(m 



a (3,3) 
2 | < , 



y(3,3) 



(logs) 



8Q!2 Tr r 



m 2 K + 



a (3,3 
16a2 



(m 2 K 



jr 



Nm 2 ss 



(61) 



In this case, 012 appears multiplied by m 2 K — m 2 (these are the tree-level 
masses, directly proportional to m s — rrid), but all higher order logarithmic 
terms proportional to a-i are subtracted in the ratio. Also, the NLO terms from 
the ©(3.3) p era tor appear in just the combinations given in Table 5. As we will 
show, the effect of the subtraction on K — > tt and K — > tttt is to eliminate the 
a 2 term, including all higher order corrections proportional to ct2 , and the NLO 
LEC's are modified to the values in Table 5, just as in the ratio for K — > 0. 
After the subtraction of K — > it, at NLO one is left with 



(n+\<D^\K+) - 2^(*+\el**>\K+) = ?pWm.(l + M - yM^roi ~ el, rot W K 

-2( e io,rot - e r S5 )m 3 K m n + (2e^ rot + 2e\ rot - 8e r 3g )m 2 K ml 

+ ( 2e 35 - e\^)ra K m\ + Ne r 2 rot m 2 K m 2 ss - Ne r 14 m K m n m 2 ss ] . (62) 

where we have indicated the coefficients that have undergone a chiral rotation to 
the form of Table 5 with a subscript, for brevity. As expected, the dependence 
on a 2 vanishes. 

The amplitudes in the partially quenched theory for K — > irir at UK1 and 
UK2 are 
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(7T + 7T-\O^\K ) ct = 8i^m 2 + 8ij W [(Ae r w + 2e r n +Ae r u -4e r 35 )m 2 

+2Nm 2 s e r u ], (63) 



for K — > 7T7T, ttik = = m (UK1), and 



x [(-2e\ + 6e r w + e\ x - Ae r 13 + Ae r 15 - 4e r 2 
"2e5 " 4e^ 5 + 8e r 39 )m K + 4e r u Nm 2 ss ] 
+l2i^m K {4L 4 -L 5 + 8L 6 + 4L s ), (64) 

for K -> 7T7T, m^i-joop) = 2m 7r(1 _ Zoop ) (UK2). 

In the partially quenched case there are, in general, additional complications 
in the calculation of the AI = 1/2, K — > tttt amplitudes due to threshold diver- 
gences leading to enhanced finite volume effects [23 OH] which require special 
care. We present the logarithmic terms for the infinite volume Minkowski space 
amplitudes at these kinematics in Appendix E. The threshold divergences are 
imaginary and vanish at NLO when the sea meson mass becomes equal to the 
pion mass (or in terms of quarks, m sea — m u = rrid) for any N > 1, see Eqs (E3, 
E5). In infinite volume Euclidean space one might expect the imaginary part 
should vanish since M Euclid — l/2(.M|i„ — M\ out ), but in lattice calculations 
the imaginary part shows up in the form of enhanced finite volume effects. In 
finite volume Euclidean correlation functions, unless the sea masses are chosen 
as stated (m sea — m u — mj), these enhanced finite volume effects are present, 
and they diverge as a power of the lattice volume. See |22l l3"%| for relevant 
calculations and discussions. 7 

As pointed out in Section 5, there exists a set of kinematics for K — > 7T7t 
where the kaon and both pions are at rest, bypassing the Maiani- Testa theorem 
on the lattice. The quark masses (m s and m u = mj) can be varied indepen- 
dently, where the weak operator inserts/removes energy to enforce 4- momentum 
conservation. We call this set of kinematics UKX, of which UK1 and UK2 are 
special cases. As pointed out by [35], there is a subtlety involved in calculating 
the UK1 kinematics, and the e prescription must be applied. One must take 
the e — > limit only after setting nix = Tn % . Given below is the most general 
expression for the LEC contribution to UKX, 



/ + -i/n(8,i)i i^0\ 4 * / , s , 8 * , 2 2 N 3m 7r (2m 7r - m K ) + ie 

(tt tt |C \K = T3ai m n (m K + m n ) + w a 2 (m K - mj ^ _ _ fe 



7 Ref 1151 states that in partially quenched lattice calculations, unless the unphysical degrees 
of freedom are above the two pion threshold, one may have serious problems with the enhanced 
finite volume effects. Again, see the note added in revision. 
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3/ 3 [4m ff (m K - m,) - it] lLV 

+ 2777^ 777^ 777^ + 2m 7r )]e£ + m| S [3(2.ZV + 8)m^ 

— 3(N + S)m,K'm 7r + (N + 6)ie]e2 — 2m K (6m n — irriKVii^ + ic)e^\ 
+ |l[( 2e io - 2e r 35 )m 3 K m^ 

+ (-2el - 4e[ 3 + 2e r 15 + 2e r 35 + ^Wk< 
+(e r n + 2e r 15 )m K ml 

+ (2e r 3 + 2e r w + e r n + 4e 13 - 4e r 35 - 8e r 39 )mi + Ne r 14 m Km7I m 2 ss }. 

(65) 

When m k = m-Ki and the limit e — > is taken in Eq (65), we recover the 
special case of UK1, given by Eq (63). One can only use UKX within the range 
m k > 777 [15] . though we show that all LEC's can still be determined. The LEC 
contribution to UKX reduces to Eq (66) at the special kinematics mx > 
and m sea = 777 u = m ( i (mss = where the imaginary threshold divergences 
vanish. 



(7r + 7r \O^ 8 -^\K ) c t = j^aim w (mK + m w ) + j^a^irnK + m^)(2m^ - ttik) 

8i 

+ p K-el + + (el -e r 5 + 2e r 10 - 2e r 35 )m 3 K m n 

+ ((-4 - N/2)el - 2el - 2e r 5 - Ae r 13 + 2e r 15 + 2e r 35 + 8e r 39 )m 2 K ml 
+ ((N/2)e r 2 + e r n + Ne r u + 2e\ h )m K ml 

+ (2e[ + (4 + N)e r 2 + 2e r 3 + 2e' 10 + e\ x + 4e 13 + Ne r u - 4e r 35 - 8e r 39 ) 



The logarithmic part of this expression is given by Eq (E10). 

The power divergent subtraction must also be performed on K — > mr ampli- 
tudes, and this requires the computation of the matrix element, (tt + tt~ |©( 3 < 3 ) \K°) . 
The subtraction to be performed is 



(n + n-\0^\K°) ee (n+n-\<D^\K°) - 2 J^(m 2 K - m%)(*+*- \Q^\K«) , 

(67) 

and the result of this subtraction at NLO is to eliminate the a.2 term, and to 
transform the NLO coefficients to the form of Table 5, just as in the case of the 
K — > 7r subtraction. This is exactly what is required, since the NLO coefficients 
that appear in Table 5 always appear in the transformed (finite) combinations 
in physical quantities, such as K — > irir at physical kinematics. 

There is a subtlety in computing K — > tttt at the kinematics where mx = 
(UK1) that must be considered when the matrix element of the O^ 3,3 ) operator 
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is computed. One expects the power divergence in (tt + tt IC^ 8 ' 1 ) |if ) to vanish 
at mx = m^ by CPS arguments, as discussed in |14|. However, a naive calcula- 
tion Of (7T + 7T" -\Q(3,3)\K°) in Minkowski space shows that a factor of mx — m v 
appears in the denominator, potentially cancelling the mx — m^ in the coeffi- 
cient multiplying (7r + 7r~ -|e( 3 ' 3 )|ir°) in Eq (67). This point was clarified in 39 . 
As they point out, it is crucial to use the e prescription in order to have a well 
defined Minkowski space amplitude. The e — > limit must be taken after the 
nix — * m-Tr limit. Thus, at leading order in ChPT, the K — > irn amplitude for 
the 6^ 3 ' 3 ' operator at UK1 is 



(7r+7r-|e^ 3 >|i^°> = lim 



4ia( 3,3 ) 3m n (2m n — nix) + * e 
k— >m„ 3/ 3 4m T (mjf — m^) — ie 



(68) 



As one can see from Eq (68) after taking the first limit, there is a pole at this 
kinematics. This pole in the denominator comes about from the graph of Fig 1, 
C2, where the kaon is annihilated by the tadpole operator. As discussed by |39| . 
it is useful to consider the corresponding amplitude in finite volume Euclidean 
space. In this case, the divergence is regulated by the finite time extent of the 
lattice, and the amplitude becomes proportional to (to LO in ChPT), which 
is the difference in time between the weak operator insertion and the two pion 
sink. When one multiplies (7r + 7r _ |0( 3 - 3 ) \K °) by m\ — m\ at exactly mx = m-*, 
then the contribution from the (~)( 3 > 3 ) term in Eq (67) is identically zero, and 
the power divergence vanishes at mx = m n , as expected from CPS symmetry. 
In the previous version of this paper, as well as in |13| . we did not properly 
appreciate this subtlety; we correct Eqs (31) and (D6) of pi] in Appendix F. 

The subtraction is necessary in UKX at all accessible values of the meson 
masses except UK1, as discussed above 8 , and at mx = 2m^ (UK2) because 
then there is no 4-momentum insertion at the weak vertex. To the extent that 
one cannot set mx exactly equal to 2771^ on the lattice, it becomes necessary to 
perform a small subtraction at this kinematics. Since the power divergences in 
UKX are proportional to mx — ^m^ [Eqs (65, E10)], the best place to investigate 
UKX numerically is in the vicinity of mx = 2m w , where we hope the power 
divergences will not be intractable. 

The NLO LEC's for the (8,1) case can be obtained as follows. From the 
subtracted K — > tt amplitude, Eq (62), one can obtain the leading order LEC, 
a\. If one uses the LEC combinations obtainable from Eq (61) for K — * 0, 
e 2,rot and e i,rot ~ e 5,rot> one can als0 obtain from Eq (62): e\ rot + e r 3 rot - 4e r 3g , 
e w rot ~ e 35J 2e^o rot — and e\ 4 . Using this information along with the 
linear combinations one can get from Eq (66) (after the subtraction) it is 
possible to obtain en + 2ei5. ro t and e\z tT0 t. Notice that UKX provides ad- 
ditional redundancy over that of UK1 and UK2 alone. For the construction 
of the physical K irir amplitude we need these seven linear combinations: 
[&2,rou &i,rot + z?>,rot — 4e3g, eio, rot — 635, 2eio, rot — en, en + 2ei5 jro t, ei3 !ro t, ei4,\. 

8 Note that UK1 is only accessible in the full theory 1151 . 
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One can verify that with these linear combinations it is possible to obtain the 
linear combinations in Eq (58). 

The logarithmic and Gasser-Leutwyler counterterm contributions to the am- 
plitudes presented in this section are given in Appendix E. 

9 Checks of the Calculations 

The logarithmic terms in the Appendixes of this paper are rather lengthy, and 
so checks are important. The first check these expressions must pass is that 
the divergences from the one-loop insertions cancel those of the divergent coun- 
terterms. This was checked for all expressions in this paper. Another check is 
that an expression reduces to some other in the appropriate limit. For example, 
in the SU(3) limit, the equations in Appendix C reduce to those of in the 
same limit, as well as those of |3(J) . modulo renormalization scheme dependent 
constants. That is, the logarithmic terms agree, but the scheme dependent m 4 
coefficients differ. 

The K — > 7T7T amplitudes in the full theory for the (8,8) 's in Appendix D agree 
with Pallante, et al. 03], as well as with ^U] (where only numerical values were 
given). Also, the PQ, K — ► 7r amplitudes of Appendix D for the (8,8)'s agree 
with JOj when they reduce to those of the full theory, in the SU(3) limit with 
"T-sea = tn va i- In the partially quenched theory, |2J has done K + — > 7r + in the 
SU(3) limit. By taking the appropriate linear combinations of the AI = 3/2 and 
1/2 amplitudes given in Appendix D we can compare to this special case, where 
we find agreement. In Appendix E, Eq (El) [PQ K ^ for the (8,l)'s] can be 
compared directly with |3L)| . where it agrees to within renormalization scheme 
dependent constants. Eq (E2) [PQ K — > it for the (8,l)'s] also agrees with [377] 
in the SU(3) limit modulo the renormalization scheme dependent constants. Eq 
(E3), K — ► 7T7T at UK1, reduces to that of the full theory for mss — m, N = 3, 
and can be compared with our previous paper ^2]- Note, however, there is an 
error in this quantity in |13| which has been corrected in Appendix F of this 
paper. Eq (E5), K — > tttt at UK2, does not reduce to that of the full theory 
since m s ^ rrid — m u , but the sea quarks were taken to be degenerate. Note 
that the logarithmic parts of the u<i term vanish for this on-shell quantity just 
as in the full theory. 

10 Conclusion 

This paper demonstrates that all of the ingredients necessary to construct all 
of the K — > 7T7T amplitudes to NLO in the full theory can be obtained without 
3-momentum insertion on the lattice, which reduces the computational cost of 
obtaining the NLO LEC's; all that is necessary in the needed K it amplitudes 
is the use of non-degenerate quark masses such that m 1 ^ ^ m^ a *. It was also 
demonstrated that all of the ingredients needed to produce e'/e to NLO are 
obtainable from partially quenched ChPT. In the case that N = 3, the LEC's 
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are those of the full theory |20I - 9 The partially quenched amplitudes were 
calculated under the assumption that both the valence and sea quark masses 
are small compared to the rf mass so that the rj' can be integrated out. This 
means that the N = limit of our amplitudes are not those of the quenched 
approximation, and this has been discussed elsewhere (3U| ■ (The terms inversely 
proportional to powers of N become quenched chiral logs.) We point out that 
we are using the PQS method [3|>|, where only the sea quarks propagate in the 
loops of Fig 3c (See Section 6.1 for the correspondence between Fig 3c and the 
traditional form of the eye-diagrams in Fig 2 for the various operators), and that 
the necessary ingredients to obtain the (8,l)'s are essentially unchanged from 
[*H?| in this prescription. 10 The PQN method, however, may not be adequate to 
determine the LEC's to NLO using only the ingredients needed of the full theory, 
except in the case where it becomes the full theory (N — 3, m sea = m va i). 

The PQChPT formulas in this paper arc valid for N = 2, however, and this 
should be useful for the work in progress by RBC with N = 2 dynamical flavors 
of domain wall quarks 0U], though the values of the LEC's determined from 
these calculations will not necessarily be those of the full theory. One would 
hope, of course, that the N dependence will not be so severe, and that this 
calculation will not be so far from the full theory. Ultimately, one would like to 
check this with a full N = 3 calculation. 

We show how the bilinear (3, 3) operator is used to eliminate the power 
divergences due to mixing with lower dimensional operators to all orders in 
(PQ)ChPT for the AI =1/2 amplitudes. The subtraction is performed to all 
orders in ChPT J5J. This is important, because the higher order power divergent 
parts can overwhelm the physical terms one is trying to calculate. 

We have pointed out that the (8,8) K — > itit amplitudes can be constructed to 
NLO using only partially quenched K — > tt amplitudes with degenerate valence 
quark masses and without momentum insertion. Also, we showed how K — > 
can be used to perform the AI = 1/2 power divergent subtraction. K — > tt 
calculations with nondegenerate valence quark masses would provide additional 
redundancy in obtaining the needed NLO LEC's. 

Finally, we point out that the threshold divergences that lead to enhanced 
finite volume corrections to the lattice calculations of AI = 1/2, K — ► tttt 
amplitudes at NLO vanish in the Minkowski space amplitudes considered in 
this paper when the sea quark mass is equal to the up and down quark masses 
{m sea = m u = rrid). This conclusion remains true in the N — 2 case for the finite 
volume Euclidean correlation functions that are relevant for lattice simulations, 
as demonstrated by [T2|- The N = 3 case had been shown to be problematic 
due to the presence of enhanced finite volume effects unless one is working in 

9 By studying finite volume Euclidean Green's functions in PQChPT, 1151 confirmed that 
one can use the partially quenched theory for AI = 1/2, K — > 7T7t amplitudes for N = 2, and 
our choice for the sea quark mass, m sea = m u ^. However, enhanced finite volume effects are 
present in the N = 3 case unless the sea quark masses are pairwise equal to the valence quark 
masses, such that one is in the full theory. 

10 The only change is that instead of UK1, one may need to use UKX, as the former has 
difficulties which were pointed out by 1151 . See also our note added in revision. 
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the full theory, with all sea quark masses equal to the corresponding valence 
quark masses ^Sl- See the note added in revision for further discussion of this 
issue. 

11 Note added in revision 

Since the original posting of version 1 of this paper on the preprint archive, Lin, 
et al. 53 have submitted a paper motivated at least in part by our version 1. 
They have done a calculation of the relevant finite volume Euclidean correlation 
functions, and they have made several important observations regarding our 
attempts to obtain K — > hit amplitudes at NLO, which we briefly review, and 
then we make some comments. They point out that the case of the AI = 
1/2 K — > 7T7T amplitudes at degenerate quark masses (what we call UK1), has 
difficulties in the full theory, and is intractable in the partially quenched theory, 
even at our special kinematics, m sea = m u> g. 

The difficulty with UK1 in the full theory is that one must disentangle various 
two meson final states in order to obtain the two pion final state. These two 
meson states have different energies, and since they appear in the correlation 
function with different exponentials in time, they will be difficult to obtain. 
See JH] for further details. They have also discovered that when unphysical 
degrees of freedom propagating in the meson re-scattering diagram are light 
enough to go on shell, they cause enhanced finite volume effects, which cannot 
be eliminated by going to larger lattices, and they make the extraction of such 
amplitudes impossible in the infinite volume limit. This problem afflicts the 
AI = 1/2 K — > irir amplitudes for degenerate quark masses (UK1) in the 
partially quenched theory, and so it cannot be used in our attempts to get the 
LEC's to NLO. 

However, it was also pointed out in |15j that the enhanced finite volume 
effects do vanish when the unphysical degrees of freedom are heavier than the 
light quark mass. This leads to the conclusion that when one is working at UKX 
(initial and final mesons at rest) if ttik is strictly greater than m v , assuming 
also m sea — TTtu.d and N = 2, the enhanced finite volume effects vanish. Thus 
UK2 (rriK = 2m,), and the kinematics points of UKX (with mx > rn w ) are 
useable in the partially quenched theory. Also, they point out that since the 
unphysical degrees of freedom cannot be lighter than the light quark mass, if 
one is working in the N = 3 theory, one must use full QCD. That is, the idea 
of varying the sea and valence quark masses independently |2(J| will not work 
for AI — 1/2, K — > 7T7T amplitudes without introducing enhanced finite volume 
effects. 

As discussed in [EJj it is not possible to use UK1 in order to obtain some 
of the NLO LEC's in the partially quenched case. We suggest that there may 
be a window where the quark masses are light enough and the lattice size is 
small enough so that the formulas of finite volume partially quenched ChPT 
can be used to extract LEC's from numerical data. Whether or not this proves 
feasible, we have found that we do not need the information from UK1 if we 
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use the K — ► irir kinematics accessible to the lattice that we call UKX. We have 
presented results at UKX (all three mesons at rest, in general, requiring energy 
insertion with mx > rn^) in section 8, of which UK1 and UK2 are special 
cases. According to there will not be enhanced finite volume effects at this 
kinematics (when N = 2, m sea = m U) d and mx > Tn„), though, in general, one 
will need to do the power divergent subtractions, as in the K — > ir case. We 
demonstrated in Section 8 that one can obtain all of the LEC's needed for the 
(8,1), K — > 7T7T amplitudes using the UKX kinematics points, along with the 
LEC's obtainable from K — > and K — > it. We conclude that it is possible 
to obtain all of the needed LEC's in the partially quenched theory for N = 2, 
though, as noted by |15j . an N — 3 determination will require the full theory 
for the information needed from lattice (8,1), K — > ttit amplitudes in order to 
avoid the enhanced finite volume effects. 
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APPENDIX A 



Appendixes B-E contain the finite logarithm and Gasser-Leutwyler counterterm 
contributions to the amplitudes presented in this paper. They were calculated 
using the FeynCalc package 01] written for the Mathematica 02] system. 
These expressions involve the regularized Veltman-Passarino basis integrals Aq, 
B and C |E|: 

Mm 2 ) = Y^yXln— , (Al) 
f 1 1 

B (q ,ml,m%)= / dx— — -[1 + ln(-x(l - x)q 2 + xm\ + (1 - x)m%) 
Jo ( 47r /j 

-ln M 2 ], (A2) 



s~i fn 2 2 2 2 2\ 

G o (0,<7 ,q ,m 1 ,m 1 ,m 2 )- 



dx- 



■x(l — x)q 2 + xm\ + (1 



x)m 



2 ' 



(A3) 
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Note that the original Veltman-Passarino integrals did not involve ChPT, 
and so the pseudoscalar decay constant, /, is not part of the original definitions 
of the integrals, but is inserted here for convenience. 



APPENDIX B 



At 1-loop order in the partially quenched theory the pseudoscalar decay con- 
stants and masses are renormalized such that J-k.k = / ^1 + — ^y- 
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i Am K \ , 
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For degenerate quark masses at 1-loop order, 'm 2 K( ^ l _ loop ^ — m 2 ^ 1 _ loop - ) 

U = fK = f(l + ^f), 
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APPENDIX C: Log Corrections to full ChPT 



The logarithmic corrections to the K — > it amplitudes in the full theory when 
3-momentum insertion vanishes are 
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These are the simplified versions of J3|> E^ 13 (C2), (D3) and (D4), respec- 
tively, when q 2 — (mK — m^) 2 . 
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APPENDIX D: PQ Log Corrections to (8,8) 's 



The logarithmic corrections for the (8,8) amplitudes relevant for the determi- 
nation of K — > 7T7T are given in this section. The logarithmic corrections to 
K — > 7T7T in the full theory were calculated first in jlOl l4*4*j , and are included 
here for completeness. 
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The A/ = 3/2, ii" — * 7r logarithmic corrections are given by 



(D2) 



(7r +| (8,8),(3/2)| X+)iog = 



4a 8 



+ 



+ 



2m K m 7r B (q 2 7 m 2 K ,ml) - N(A (m 2 sS ) + 3A {ml s )) 
2mV 



m* + m 2 K (m 2 



SS 



N(2mj 



im 2 K m 2 + ml) 



A (m 33 ) 



2ml) 



K mii-K 

mj + m 2 K (m 2 s 

Nml(ml — m 2 K ) 
A/ K A// 



M™1) + 



(m 



/ 



l ss) 



(D3) 



For the AI — 1/2, K — > 7r corrections there are (at least) two possibilities, 
when the electroweak operator is partially quenched and when it is not. The 
following amplitude corresponds to quenching the short distance electroweak 
operator, neglecting the type of contraction in Fig 3c (with a photon or Z 
replacing the gluon), 
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while the next corresponds to where the short distance electroweak operator is 
not quenched, and valence quarks do propagate in the loops of Fig 3c (again 
with a photon or Z replacing the gluon) . 
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In the case of degenerate valence quarks the above expressions reduce to 
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Note that the two AI =1/2 expressions reduce to the same thing in the SU(3) 
limit. 



APPENDIX E: PQ Log Corrections to (8,1) 's 

The logarithmic corrections for the quantities relevant for the determination of 
the (8,1), K — > 7T7T amplitudes are given in this Appendix. The logarithmic 
corrections to the physical K — > irn amplitude have been done by |2HH2H], and 
we refer to ^3]> Eq. D10, for the amplitude in our conventions. 
The logarithmic corrections to K — > and K — > it are given by 
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In the case of degenerate valence quarks, the above expression becomes 
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The K — > 7T7T amplitude at UK1 in infinite volume Minkowski space is given 



by 
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for nix = = m, where 
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Expression E5 is real for m 2 ss > m 2 . When m 2 ss < m 2 , E5 has an imaginary 
part. The K — > tttt amplitude at UK2 in infinite volume Minkowski space is 
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for m K (i-i oop ) = 2m 7r (i-i oop ), where 
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Note the imaginary threshold divergences in both (E4) and (E6). On the 
lattice they are expected to contribute in the form of enhanced finite volume 
effects. See, for example, When mss < mx/2 = m n , Eqs. (E7) and (E8) 
have imaginary parts. 

Eq (E6) is most useful in fits to lattice data at the special kinematics m sea — 
m u,d i m ss = m 7r)j N = 2. In this case it reduces to 
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The logarithmic contribution to UKX (kaon, pions at rest) in the special 
case of m sea = m u = m^ (mss — m n) & n d mx > m^ is given below. In this 
expression, as in all others in this set of Appendixes, q 2 — (mjf — m v ) 2 . 
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APPENDIX F 



The absence of the a 2 terms in K — > mr at UK1 requires some corrections to 
, presented here. Eq (31) of 13 should be 



(n + n-\O^\K ) ct = 8z^m 2 +8z^[44 + 2e^+44 5 -4eS 5 ]. 
In Appendix D, Eq (D6) of jT^], the correct equation should read 
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